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Abstract 

A natural connection between rational functions of several real or complex vari¬ 
ables, and subspace collections is explored. A new class of function, superfunctions, 
are introduced which are the counterpart to functions at the level of subspace col¬ 
lections. Operations on subspace collections are found to correspond to various 
operations on rational functions, such as addition, multiplication and substitution. 
It is established that every rational matrix valued function which is homogeneous 
of degree 1 can be generated from an appropriate, but not necessarily unique, sub¬ 
space collection: the mapping from subspace collections to rational functions is 
onto, but not one to one. For some applications superfunctions may be more im¬ 
portant than functions, as they incorporate more information about the physical 
problem, yet can be manipulated in much the same way as functions. Previously 
subspace collections had been introduced when there was an inner product on the 
vector (or Hilbert) space, and appropriate subspaces were mutually orthogonal. In 
that setting certain normalization and reduction operations on subspace collections 
led to a continued fraction expansion of the associated function, which allowed one 
to bound the function in terms of a set of weight matrices and normalization ma¬ 
trices that are derived from series expansions. Here we also initiate the theory of 
normalization and reduction operations, appropriate when there is no inner product 
on the space. 


1 Introduction 

This Chapter 7 of the book ’’Extending the Theory of Composites to other Areas of Sci¬ 
ence” , edited by Graeme W. Milton, is concerned with developing the theory of subspace 
collections, particularly nonorthogonal subspace collections. Subspace collections have a 
rich algebraic structure, and a close connection with rational functions of several real or 
complex variables. Here we are interested in three types of subspace collections: first, 
finite dimensional vector spaces "H that have the decomposition 

'hi = Li (B £ ® = 'Pi © P 2 © • • • © Pn.) (1-1) 
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which we call a Z (n) subspace collection; second hnite dimensional vector spaces /C (over 
the real or complex numbers) that have the decomposition 


/c = ^©dr = v©Pi©p2 




( 1 . 2 ) 


which we call a Y(n) subspace collection, where the S and entering (1.2) are not to be 
confused with the subspaces £ and J' entering (1.1); and third hnite dimensional vector 
spaces /C (over the real or complex numbers) that have the decomposition 

JT = © Pi © p2 © • • • © 


IC = £ 


(1.3) 


which we call a superfunction F^{n). In a superfunction the space and the space 
are called the input and output subspaces respectively, and they have the same 
dimension. For superfunctions we require the technical condition that for any choice of 
vectors E^, G and G there exist vectors E G and J G JT” such that 

E^ = n^E, = n^j, E° = n^E, = n^j. (1.4) 


As we will see there is a very close direct connection between a superfunction F^{n) 
and a Y (n) subspace collections, and also many connections between them and Z[n) 
subspace collections. All are intertwined and that is the beauty of the theory. Z{3) and 
Y{2) subspace collections, and superfunctions P®(1) can be visualized in 3-dimensional 
space, and examples of these are given Figure [Tj 



Figure 1: Shown in (a) is an example of a Z{3) subspace collection, in (b) a Y (2) subspace 
collection, and in (c) a superfunction P^(l). The rays denote one-dimensional subspaces: 
they should really be drawn as lines, but for clarity they are drawn as rays and should 
be extended in the opposite direction as the ray. The circles, which look like ellipses as 
they are tilted, represent two-dimensional subspaces. 
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One reason V(n) snbspace collections, Z(n) subspaces collections, and superfunctions 
F*(n) are important is because they arise in many physical problems. For examples in 
network theory and in the theory of the effective moduli of composite materials, see 
the review in Chapter 2 of this book (Milton 2016) and Milton (2002). There are also 
many other physical problems where subspace collections arise as is apparent in Chapters 
1,3,8,9, 12, 13, and 14 of this book (Milton 2016). In physics applications the subspaces 
are usually orthogonal with respect to some inner product on the space V. or JC but as 
this chapter shows the theory of them can be developed without reference to an inner 
product. This generalization is important to make contact between general rational 
functions of complex variables, thus extending the notion of a function: hence the name 
superfunction. The generalization is also important for applications, such as speeding 
up numerical methods for calculating the helds that solve the problem: we will see an 
example of this in the next chapter. 



Figure 2: Two routes to solving a physical problem formulated in terms of subspace col¬ 
lections. It is suggested that the route on the right may result in a better approximation 
as more information is kept. 

It may very well be the case that superfunctions become more important than func¬ 
tions in some applications, as suggested by the flow chart of Figure The reason is 
that when one extracts the function from a superfunction, which we will see how to do 
shortly, one generally loses information that is contained in a superfunction. For exam¬ 
ple, in the context of physical problems where there is an inner product on the space 
this information may came in the form of a series expansion for the helds up to a given 
order, and from this series expansion one can extract the “weight matrices” and “normal- 
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ization matrices”, introduced by Milton and Golden (198^ and Milton (1987a, 1987b). 
These matrices basically encode the information about the “angles” between the various 
subspaces (when there is an inner product). One can then develop a continued fraction 
expansion for the function associated with the superfunction, with the normalization 
factors and weight matrices that enter it at each level having the property that they are 
positive semidehnite, with the weight matrices summing to one. Truncating the contin¬ 
ued fraction gives approximations to the function, that are similar in some respects the 
diagonal Fade approximants, and in fact give bounds on the function if the truncation is 
done appropriately. The information contained in the weight matrices and normalization 
matrices, cannot in general be recovered (at least when n > 4) from the series expansion 
of the associated function. (Although one can potentially determine these matrices from 
the series expansion of the functions associated with coupled held problems, as shown 


in Chapter 9 of this book (Milton 2016|). This theory was established by Milton and 


Golden (1985) and Milton (1987a, 1987b[ [1991 ). 


Milton (2002)) for the case of Z{n) 


L6|). 


see 


also Chapters 19, 20 and 29 in 
subspace collections, for any integer n > 1. In this 
paper develop the basic theory of subspace collections in the case where there is no inner 
product on the vector space 77 or /C. We also make the hrst steps towards generating 
continued fraction expansions in the case where there is no inner product on the vector 
space 77 or /C. 

Let us hrst suppose V and U are one-dimensional. We will see that there are generally 
homogeneous (of degree 1) rational functions Y (zi, ^ 2 , • • •, Zn) and Z{zi, 2 : 2 , • • •, Zn) (over 
the real or complex numbers) of degree 1 that are associated respectively with these 
Y{n) and Z{n) subspace collections, where Z{zi,Z 2 , ■ ■ ■, Zn) satishes the additional con¬ 
straint that Z(l, 1,..., 1) = 1. Conversely, we will see that given any rational functions 
Y(zi, Z 2 , ■.., Zn) and Z{zi, Z 2 , ■ ■ ■, Zn) with these properties, then there exists at least one 
subspace collection realizing these functions as its associated function. There are also 
operations on these subspace collections that correspond to operations on the associated 
function, such as substitution. 

For superfunctions the simplest case is when the input and output spaces and 
are one-dimensional. Then with a specihc basis for and the corresponding func¬ 
tion F(zi, Z 2 , ■ ■ ■, Zn) is 2 by 2 matrix valued with the elements Fii{zi, Z 2 , ■ ■ ■, Zn) and 
7 ^ 22 ( 2 ^ 1 , Z 2 ,, Zn) being homogeneous of degree zero, the element i^i 2 (^i, ^ 2 , • • •, Zn) be¬ 
ing homogeneous of degree minus 1, and 7 ^ 21 ( 2 ^ 1 , ^ 2 , • • •, Zn) being homogeneous of degree 
1. There are operations on superfunctions that correspond to addition, multiplication and 
forming an inverse (and hence division) of the associated functions. So superfunctions 
form an algebra. Also one can do substitutions at the level of subspace collections. Ac¬ 
tually the operation of addition of superfunctions are naturally done with the associated 
Wproblem, although one could equally do them with the associated inverse F-problem 
(where the spaces T and J7 are interchanged). Thus there is an inherent ambiguity of how 
one wants to dehne addition of superfunctions. The dehnitions of addition, multiplica¬ 
tion. and substitution of subspace collections may seem a little complicated and abstract, 
yet they are the exact counterpart of similar operations one may do on multiterminal 
electrical networks, and they do produce the corresponding action on the associated func¬ 
tions. (In fact it was thinking about electrical circuits which guided the construction of 
these operations in a more general setting). 

When V and 77 have dimension greater than 1, then Y(zi, Z 2 ,..., Zn) and Z(zi, Z 2 ,..., Zn) 
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get replaced by linear operator valued functions 2 : 2 , • • •, Zn) and Z( 2 :i, 2 : 2 , • • •, Zn) 

which map V to V and U to U respectively. Similarly, the function F(2:i, 2 : 2 ,..., Zn) 
should really be thought of as a linear operator mapping to 

The original motivation for studying subspace collections, and their associated func¬ 
tions, arose from the study of the effective conductivity tensor Z of periodic com¬ 
posite materials. For a composite with n isotropic phases, with scalar conductivities 
2 : 1 , 2 : 2 ,. ■ ■, Zn, the effective conductivity tensor was found to be a homogeneous (of de¬ 
gree 1) analytic function Z(zi, Z 2 , ■ ■ ■, Zn) of the component conductivities with positive 
dehnite imaginary part when the component conductivities have positive imaginary part 
Bergman 1978 Milton 1979, 1981a, Golden and Papanicolaou 1983| (see also Chapter 


18 of Milton (2002)). It was also recognized (Milton 1987a, 1990) that the problem of 
determining the effective conductivity function could be formulated in terms of three mu¬ 
tually orthogonal spaces in the Hilbert space "H of square integrable functions: namely 
the space lA of constant helds, the space £ of periodic square integrable electric helds 
(having zero curl), and the space J of square integrable current helds (having zero di¬ 
vergence), and if the composite had n isotropic phases, with conductivities 2 : 1 , 2 : 2 ,..., Zn, 
then it was also natural to decompose 77 into the direct sum of n mutually orthogonal 
subspaces Vi,V 2 , ■ ■ ■ ,Vn where Vi consists of those square integrable helds which are 
nonzero only within component i. This formulation, in terms of a Z[n) subspace col- 


lection, evolved out of earlier Hilbert space formulations of the problem ( 

Fokin 1982 

Kohler and Papanicolaou 19821 Papanicolaou and Varadhan 19821 

Golden and Papan- 

icolaou 1983[ 

Kantor and Bergman 1984| 

Dell’Antonio, Figari, and Orlandi 1986 

) and 


can easily be extended to the elastic, thermoelastic, piezoelectric, and poroelastic equa¬ 


tions of multiphase and polycrystalline materials (see, for example. Chapter 12 in Milton 


(2002)). The formulation has proved to be particularly important in the theory of exact 


relations of composite materials ( 

Grabovsky 1998 

Grabovsky and Sage 1998 

Grabovsky 

land Milton 1998 Grabovsky, Milton, and Sage 2000; 

Grabovsky 2004) (see also Ghapter 


17 in Milton (2002)) where one seeks microstructure independent relations satished by 
ehective tensors. For two-dimensional polycrystals a complete correspondence was estab¬ 
lished between subspace collections and a representative class of multiple rank laminate 
polycrystal geometries (Clark and Milton 1994), thus showing that the subspace collec¬ 
tion of any two-dimensional polycrystal, with any conhguration of crystal grains, could 
be approximated arbitrarily closely by the subspace collection of one of these multiple 
rank laminate polycrystal geometries. 

Curiously the connection between Z{n) subspace collections and the ehective conduc¬ 
tivity allowed the ehective conductivity function Z{zi, Z 2 , ■ ■ ■, Zn) to be expanded as a new 
type of continued fraction, involving matrices of increasing dimension as one proceeds 
down the continued fraction when n > 2 (Milton 1987a, 1987b, 1991 see also Chapters 
19, 20 and 29 in Milton 2002). The coefficients in the weight and normalization matri¬ 
ces entering the continued fraction can be expressed in terms of inner products between 
helds that enter the series expansion of the solution held in a nearly homogeneous medium 
(with all the conductivities 2 : 1 , 2 : 2 ,..., being close to one another). One application of 
the continued fraction expansion has been to obtain bounds on the diagonal elements of 
the complex ehective conductivity tensor of a three phase conducting composite, with 
complex conductivities Zi, Z 2 and z^, that were tighter than bounds obtained by any 
other method (see hgure 4 in Milton 1987b). This procedure essentially extended to 
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multivariate functions the procedure, using successive fractional linear transformations, 
that was used to obtain bounds ( |Baker, Jr. 1969 ) on the values in the complex plane 
that Stieltjes functions! can take when a hnite number of Taylor series coefficients are 
known (see also Golden and Papanicolaou 1983; iBergman 1986) where essentially the 


same transformation is used to derive bounds on the complex dielectric constant of two 


component media using series expansion coefficients, as noted in the appendix in Milton 


1986, and see [Milton 1981b where the same set of bounds is derived using a different 


procedure, namely the method of variation of poles and zeros.) 

In the case n = 2 the continued fraction reduces to a usual continued fraction expan¬ 
sion, like those continued fractions associated with Fade approximants (see Chapter 4 of 
Part I of Baker, Jr. and Graves-Morris 1981). Y{n) subspace collections enter, for exam¬ 
ple, if one eliminates from the Hilbert space the constant helds and then reformulates the 
conductivity equations in terms of the remaining helds: the driving helds are then helds 


which are constant in each phase, but have zero average value (see Chapter 19 in Milton 
2002 and references therein). The interrelationship between Z{n) subspace collections 


and Y(n) subspace collections is what gives rise to these novel continued fractions. 

Finite dimensional Z{n) and Y{n) subspace collections also arise naturally in the 
study of the ehective resistance of electrical circuits constructed from n types of resistors 
having conductances Zi, Z 2 , ■ ■ ■ Zn (see Chapter 20 in Milton 2002). This is not surprising 
as periodic resistor networks can be seen as discrete approximations to conducting com¬ 
posite materials (see, for example, Milton 1981a and Figure 8.5(a) in this book Milton 


2016[). Figure 1^ illustrates a discrete network of impedances, and gives an indication of 


the physical meaning of the Z{n) and Y{n) subspace collections in this context. 

In this hgure, the vector space H is 6 -dimensional, and is the direct sum of the two- 
dimensional space Vi consisting of helds that are nonzero only along the resistors CiZi 
and C 3 Z 1 ] the two-dimensional space V 2 consisting of helds that are nonzero only along 
the resistors C 2 Z 2 and C 5 Z 2 ; and the one-dimensional space V 3 consisting of helds that are 
nonzero only along the resistor 04 ^ 3 . The response of the network, when one terminal is 
grounded (with zero voltage) is a 3 x 3 matrix. When it acts on the vector, having as 
elements the voltages at the three remaining terminals, it gives the three currents howing 
through these terminals. The 3x3 matrix valued function Z(zi, Z 2 , Z 3 ) gives the matrix 
valued response relative to the response when zi = Z 2 = Z 3 = 1. Now, let us imagine all 
the resistors, or impedances, in (a) are on one side of the circuit board, with the terminals 
being conducting posts that penetrate the board. On the other side of the board these 
posts are connected to a tree-like graph of batteries (or alternating current sources if 
the helds vary sinusoidally in time) shown in (b). The three helds in these batteries 
constitute the space V. The Y (3) subspace collection contains helds on both sides of 
the board, in /C = H © V. The associated 3x3 matrix valued F-function Y{zi,Z 2 , Z 3 ) 
gives the current going through the three batteries, in response to the voltages across 
them. Note that Y {zi,Z 2 , Z 3 ) is not diagonal: a voltage across one battery, sends current 
through the other two batteries, even when they have zero voltage across them. 

Superfunctions are a natural generalization of multiport electrical circuits with input 
ports and output ports, as illustrated in Figure]^ The function F gives the currents and 
potential drops across the output batteries/resistors that are generated in response to 
currents and potential drops across the input batteries. Note that the networks associated 
with superfunctions automatically satisfy the “port condition” that the net flow of current 
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Figure 3: Shown in (a) is a 4 terminal electrical network, which is representative of a Z{3) 
subspace collection. Here the Cj are real positive scaling constants: the conductance of 
each element is CjZk where Zk is real or complex (when is complex we should refer to cjZk 
as an admittance rather than as a conductance). Complex values of 2 ; are appropriate 
when the applied potentials vary sinusoidally with time, and some of the impedence 
elements are capacitors or inductors. Figure (b) shows the batteries on the back side of 
the circuit board, representing the space V, which combined with the resistors on the 
front side is representative of a 4^(3) subspace collection. The F-function Y{zi,Z 2 ,z^) 
gives the current going through the three batteries, in response to the voltages across 
them. 


from the input terminals to the output terminals is zero. 

In this chapter we show that the connection between hnite dimensional Z{n) and 
Y(n) subspace collections and homogeneous (degree 1) operator valued rational functions 
Y(^ 1 , Z 2 ,..., Zn) and Z(zi, Z 2 , ■ ■ ■, Zn) persists even when the subspaces in each decom¬ 
position are not necessarily mutually orthogonal, and indeed even in the absence of an 


is no inner product. Accordingly some steps in the analysis, and some assumptions, need 
to be revised. In this more general setting we can generate, from an appropriate Z(n) 
subspace collection, any desired scalar valued rational function Z(zi, Z 2 , ■ ■ ■, Zn) satisfying 
the homogeneity property Z(l,l,...,l) = 1. 

It is to be emphasized that subspace collections, with the associated rules for addition, 
multiplication and substitution, are algebraic objects in their own right: there is no need 
to think of the associated analytic functions (that are in general operator valued), except 
that the correspondence makes it easier to think about subspace collections. The resistor 
network examples of Y (n) subspace collections made it possible for me to see how the 
operations of addition, multiplication and substitution of subspace collections should be 
dehned in the general case. 

My belief is that the geometrical structure of subspace collections (and in particular 
superfunctions) will be reflected in the algebraic geometrical structure of their associated 
rational functions. If this is the case, understanding the topological features of subspace 
collections might shed light on the geometrical features of algebraic varieties. While this 
paper does not directly address this issue, it sheds the hrst light on the relation between 
hnite dimensional subspace collections and rational functions of several complex vari¬ 
ables, in the case where the subspaces are not mutually orthogonal, and it introduces 


inner product (on the space "H or /C). The results developed in (Milton, 1987a, 1987b 


1991 and in Chapters 19, 20 and 29 of Milton, 2002) are extended to the case where there 


7 


















Figure 4: Shown in (a) is a 5 terminal electrical network, which is representative of a 
Z(3) subspace collection. Here the q are real positive scaling constants: the admittance 
of each element is CjZk where Zk is real or complex. Figure (b) shows the batteries on the 
back side of the circuit board, representing the space V, which is divided into the input 
space V^, consisting of those vectors in /C that are nonzero only in the batteries Ji and 
I 2 and the output space V^, consisting of those vectors in /C that are nonzero only in the 
batteries/resistors 0i and O 2 . Figure (c) shows a 6 terminal electrical network, and the 
naturally associated subspace V represented by the batteries in Figure (d). To convert 
this to a problem where the dimension of V is even we remove the battery at the top, 
and accordingly reduce the dimension of both V and by one. Figure (e) shows the 
input space , consisting of those vectors that are no n zero only in the batteries Ji and 
I 2 and the output space V^, consisting of those vectors in /C that are no n zero only in the 
batteries/resistors 0i and O 2 . 


superfunctions. The functions derived from superfunctions are well studied and have 
widespread applications in signal processing, control theory, network synthesis and de¬ 
sign, and in optics, acoustics and elastodynamics (usually in layered media), where they 
are called a variety of names including transfer matrices, transmission matrices, transfer 
functions, system functions, and network functions. In these contexts it is the function 
that is studied, but people do not think of the superfunction. I thank Aaron Welters and 
Mihai Putinar for drawing my attention to the connection between transfer functions 
and response functions (such the effective conductivity tensor of composites). 

We remark that for Z (3) orthogonal subspace collections, with U being one-dimensional, 
it is still an open and intriguing question as to whether there could be a one-to-one corre¬ 
spondence between them (assuming they are pruned as described in Section 15 and mod¬ 
ulo trivial equivalences between subspace collections) and scalar functions Z{zi, Z 2 , Z 3 ) 
satisfying the homogeneity, Herglotz and normalization properties. The Z-problem de¬ 
scribed the next section provides a nonlinear map from the Z{3) orthogonal subspace 
collection to an associated scalar function Z{zi, Z 2 , Z 3 ) satisfying the homogeneity, Her¬ 
glotz and normalization properties, but the question is whether one can uniquely re¬ 
cover the pruned subspace collection, modulo trivial equivalences, given only the function 







Z(zi, Z 2 , Zs)? The intriguing counting argument given in Section 29.2 of Milton (2002) 
suggests the possibility of a one-to-one correspondence. There is a similar counting ar¬ 
gument for nonorthogonal subspace collections given in Section 18, but in this case we 
will see in an explicit example that a one-to-one correspondence does not hold. 


2 Subspace collections and their associated functions 


Let /C be a vector space which has a decomposition into two different direct sums of 


subspaces 

/c = £ej = ven, ( 2 . 1 ) 

where "H itself is a direct sum of n subspaces 

l-L = Vi®V2®---®Vn. (2.2) 

Any vector K G /C has a unique decomposition into component vectors, 

K = E + J = v + H, H = Pi -fP2 + ■■■ + ?„, (2.3) 

each in the associated subspaces: 

E G J G J”, V G V, H G "H, Pj G "Pi for i = 1,2, ..., n. (2.4) 

This decomposition serves to dehne projection operators Pi and r 2 onto £ and J, projec¬ 
tion operators IIi and 112 onto V and "H, and projection operators Aj onto the subspaces 
Vi- By dehnition we have 

E = PiK, J = r2K, V = HiK, H = HsK, P^ = A,K. (2.5) 


Associated with this subspace collection is an linear operator valued function Y(zi, Z 2 , ■ ■ ■, Zn) 
acting on the space V, which is a homogeneous function of degree 1 of the n complex 
variables ^i,^ 2 , ■ ■ ■ ,Zn- To obtain the function we take each held Ei G V and look for 
vectors J and E that solve the equations 

E G J G , J 2 = LE 2 , where J 2 = n 2 J, E 2 = II 2 E, (2-6) 


with El = IIiE, where 

n 

L = (2.7) 

i=l 

We call this problem the E-problem. The associated operator Y, by dehnition, governs 
the linear relation 

Ji = —YEi, where Ji = IIiJ. (2.8) 

A necessary condition for Ji to be uniquely dehned given Ei is that 


vnj = 0, 


(2.9) 


since if J and E solve (2.6) so too will J -|- v and E, for any v G V ft 77. The inverse 
Y-problem is to solve (2.6) for each held Ji = IIiJ G V. A necessary condition for Ei to 
be uniquely dehned given Ji is that 


vn£ = 0. 


( 2 . 10 ) 
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If vi, V 2 ,..., Vm is a basis of V, then the operator Y can be represented by a matrix, the 
Y-matrix, also denoted by Y with elements such that 


Yvfc = ^ (2.11) 

i=l 

If m is even and V has the decomposition 

V = v^©v^, (2.12) 

where and have the same dimension (m/2) then we have a superfunction F^. The 
superfunction is the collection of subspaces and there is a function F associated with it. 
The helds Ei and Ji have the unique decomposition 

Ei = E^ + E®, Ji = J^ + J°, (2.13) 

with 

E^J^gV^ E°,J^eV°, (2.14) 

where the superscripts I and O refer to input and output respectively. We write 


E^ = n^Ei, fP = n^Ei, y = n^Ji, j® = n^Ji, (2.15) 


which dehnes the projections 11^ and 11^ onto the input and output spaces. Now the 
relation (2.8) can be written as 



/Yii (F^\ 

Y ®°) l^E ®) ’ 


and manipulated into the form 



( 2 . 16 ) 


( 2 . 17 ) 


which dehnes the linear operator valued function 


F = 



-{Y^OyiYii 

YOO/yioy^Y^^ — Y^^l 


_(Y/0)-i 

yoo 


( 2 . 18 ) 


which, provided the operator Y^^ is nonsingular, is the function associated with the 
superfunction. This relation can be inverted to yield Y in terms of F, 

( (F^-^)-iF^^ -(F^^)-i \ 

^ I [F'^‘^(F^'^)“^F®'® — F*^^] _F'^'^(F®‘^)“^y ’ 


provided the operator F®*^ can be inverted. The superfunction problem is for given input 
helds E^ and to hnd helds E and J that solve the Y-problem (2.6) and (2.7), with 
n^E = E^ and II^J = J^. It may happen that the superfunction problem has a solution 
when the Y-problem does not (this happens when F'®'^ is singular), and conversely the 
Y-problem may have a solution when the superfunction problem does not (this happens 
when Y^'^ is singular). 
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Another association between subspace collections and functions conies if a vector 
space l-i has the decomposition 


— lA ® £■ ® J' — T’l © 1^2 © • • • © 'Pni 


( 2 . 20 ) 


where £ and J are not to be confused with the spaces in (2.1). Any vector H G "H has 
a unique decomposition into component vectors, 

H = U + E + J = Pi+P2 + ''‘ + Pn) (2-21) 

each in the associated subspaces: 

u G W, E G J G JT”, Pj G Pi for i = 1, 2,..., n. (2.22) 


This decomposition serves to define projection operators Pq, Pi and P 2 onto W, £ and 
J , and projection operators A* onto the subspaces Vi. Associated with this subspace 
collection is an linear operator valued function Z(2;i, ^ 2 , • • •, Zn) acting on the space U, 
which is a homogeneous function of degree 1 of the n complex variables zi, Z 2 , ■ ■ ■, Zn- To 
obtain the function we take each vector e G W and look for vectors j, J and E that solve 
the equations 

n 

JgW, EgT, JgJT”, j + J = L(e + E), where L = ZiXi. (2.23) 

i=l 

We call this problem the Z-problem. The associated operator Z, by definition, governs 
the linear relation 

j = Ze. (2.24) 

If ui, U 2 ,..., Um is a basis of W, then the operator Z can be represented by a matrix, also 
denoted by Z with elements Zi^ such that 


Zu/j ^ ^ 


(2.25) 


Z=1 


When zi = Z 2 = ■ ■ ■ = Zn = ^ (2.23) has the trivial solution 


j = e, J = E = 0, 


(2.26) 


and so we deduce that 


Z(1,1,...,1)=P (2.27) 

The inverse Z-problem is to solve the equations (|2.23) for each given vector j G 7/. 


3 Some simple examples 

Consider a Y (n) subspace collection 

K, = £®£f = Y® V\ © V 2 © • • • © Vny (3-1) 
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where S,V,Vi,V 2 , ■ ■ - Vn are all one-dimensional, and is n-dimensional. Choose, as 
our basis for 1C, n + 1 vectors po G V, and p* G Vi, i = 1, 2,... n. Vectors E G and 
J G JT” can be expanded in this basis: 


E = '^E.pi, 
i=0 


n 

j = 'Y^JiPi. 

i=0 


(3.2) 


The relation II 2 J = Ln 2 E implies 


Ji = ZiEi. (3.3) 

Let us suppose that Eq = 1. Then Ei and E 2 are determined by the orientation of the 
one-dimensional subspace £ with respect to the subspaces V, Vi, V 2 , ■ ■ ■ Vn- Also since J 
has codimension 1 , there exist constants ITo, Wi,... Wn, determined by the orientation 
of the n-dimensional subspace J with respect to the subspaces V,Vi,V 2 , ■ ■ - Vn such that 

n 

J2W^J^ = 0. (3.4) 

i=0 

Let us suppose that ITo = 1. Then we have 

n n 

Jo = - ^ WiJi = - ^ WiEiZi, (3.5) 

i=l i=l 

which since Eq = 1 implies Jo = —YEq, with 

n 

Y = '^^aiZi, where ctj = HJEj. (3.6) 

i=l 

As the Ei and Wi are arbitrary constants, we see that Y can be any linear combination 
of the Zi- In particular, with WiEi = 1 and WiEi = 0 when i 7 ^ 1 we obtain 

V = ;^i. (3.7) 

As a second example consider a V(l) subspace collection 

]C = £®J = V®Vi, (3.8) 

where all the spaces £, J, V, and V\ are all two-dimensional. Choose as our basis 

for /C two vectors pi and p 2 in V and two vectors ps and p 4 in V\. Then since £ is 

two-dimensional, there generically exist constants 613 , 614,623 and 624 such that 

Pi + 613 P 3 -f 614 P 4 G T, P 2 + 623 P 3 + 624 P 4 G £. (3-9) 

Also since J is two-dimensional, there generically exist constants ^ 31 , ^ 32 , j 4 i and ^42 such 
that 

Ps + J 31 P 1 + J 32 P 2 G 77 , P4 + J 41 P 1 + e42P2 G J. (3.10) 
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So the y-problem is solved with vectors 


E = Pi + ei3P3 + ei4P4, 

El = Pi, E 2 = ei3P3 + ei4P4, 

J 2 = 2:i(ei3P3 + ei4P4), 

j = 2:i[ei3(P3 + jsiPl + j32P2) + ei4(p4 + j4lPl + e42P2)], 

Jl = ^l[(ei3j31 + ei4j4l)pi + (ei3j32 + e42j42)p25 (3.11) 

and is also solved with vectors 

E = P 2 + e23P3 + e24P4, 

El = P 2 , E 2 = e23P3 + e24P4, 

J 2 = 2:i(e23P3 + e24P4), 

J = 2:i[e23(P3 + jsiPl + j32P2) + e24(P4 + j4lPl + e42P2)], 

Jl = [(623^31 + e24j4l)Pl + {^21,jz2 + e24j42)P2- 

From these equations in follows that Y{zi) in this basis is the 2 by 2 matrix 

Y( 2 ;i) = ziA, with A = , 

V®21 «22 J 

where 


(3.12) 

(3.13) 


dll — 613^31 + ei4j4i, ai2 — ei3j32 + 642^42, 

021 = 623^31 + 624 ^ 41 , 022 = 623^32 + 624^42 • ( 3 - 14 ) 


As the coefficients 613 , 614 , 623 , 624 ,^ 31 ,^ 32,^41 and ^42 can be any complex numbers we 
desire it follows that we can realize any desired complex matrix A. By taking to be 
the one-dimensional space spanned by pi and taking to be the one-dimensional space 
spanned by p 2 we obtain a superfunction where the associated function takes the 
form 


F(;^i) 


f hii hi2lzi\ 
\&212 :i &22 / ’ 


(3.15) 


in which the parameters 611 , 612 , &21 and 622 can be any complex numbers we choose. 
As a third example consider a Z{2) subspace collection 


n =u ®s ® J = 'Pi®V2, 


(3.16) 


where the subspaces U, S, J and V 2 are all one-dimensional, while "Pi is two-dimensional. 
Choose, as our basis for "H, 3 vectors Uq G W, Eq G and Jo G 57”, and take a vector P 
as a basis for V 2 - The coefficients P^/, Pe and Pj in the expansion 

P = Pi/Uo-l-P eEq- l-PjJo (3-17) 

determine the orientation of P 2 with respect to the subspaces 7/, E- and J. In the basis 
Uo, Eq, and Jo the equations 


e -|- E — Q -|- ciP, j + J — ^iQ + Z 20 tP^ 


( 3 . 18 ) 
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with 


take the form 


e, j G W, E G J G 


Q ^ ”^1) 




~h 



? 



+ 



? 


( 3 . 19 ) 


( 3 . 20 ) 


and since Q G "Pi there exist constants Wjj^ We and PEj, which determine the orientation 
of Vi with respect to W, £ and J, snch that 


WuQu + WeQe + WjQj = 0. 
Hence we obtain the eqnations 


(3.21) 


Wue + WeE = aiWuPu + WePe + WjPj) = aW ■ P, 

Q = zi{E - aPE) + Z2aPE, 

j = Zi{e — aPu) + Z20iPu- (3.22) 


Eliminating E and a from these eqnations gives j = Ze, with 

^ _ _ (^2 — ^l)PEi7Pj7 _ 

W ■P + WEPEiz2-Zi)/z,' 

In particnlar if the snbspaces are oriented so that 


(3.23) 


W ■ P = WePe = -WuPu. 


(3.24) 


then (3.23) gives 


Z = zllz 2 , 


which with Z 2 = \ prodnces the fnnction z\ and with zi 
Also, with WePe = 0 we obtain 


(3.25) 


1 prodnces the fnnction 1 / 2 : 2 . 


Z 


Zl + 


{Z 2 — Zi)WuPu 

W P 


(3.26) 


which is a “weighted average” of Zi and 2 : 2 , Z = WiZi + W 2 Z 2 with “weights” Wi and W 2 
that snm to 1 bnt which are not necessarily positive, nor even real. 


4 Formulas for the associated functions 

Following Section 12.8 of Milton (2002) a formnla for the effective tensor Z resnlts by 
applying the operator Pq + r 2 (which projects on the space W © J7) to both sides of the 
constitntive law e + E = L^^(j + J). Solving the resnlting eqnation. 


e = 


= (ro + r2)L-i(ro + r2)(j + J) 


(4.1) 
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for j + J gives 


j + J = [{To + T2 )L-\To + r 2 )]-'e, (4.2) 

where the last inverse is to be taken on the snbspace U ® J- By applying Fq to both 
sides of this eqnation we see that 


z = r„|(r„ + r2)L-‘(r„ + r2)|-'r„, 


( 4 . 3 ) 


which is the resnlt given in (12.59) of Milton (2002). 

Another formula for Z follows from noting that for any arbitrary constant Zq ^ 0, 


[zqI — ri(L — zoI)](6 + E) — zqG + zqE — FiJ — zoFiE — zqG. 

Solving this for e + E gives 

e + E = zqIzqI — ri(L — 2;oI)] 
and applying FoL to both sides yields 

j = ^oroL[zoI ~ Fi(L — ^qI)] ^6- 
Thus we have a formula for the Z operator, 

= zqTqL[zq\ — Ti(\j — zq\)\ ^Fq = ^qFo + 2:oFo(L — zoI)[2:oI ~ Fi(L — zqI)] ^Fq, 
where we have used the identity 

To = ~ ri(L — 2;oI)] ^Fq, 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


obtained by applying Fq to both sides of (4.5). This formula (4.7) is a special case of the 


formula (12.60) given in Milton (2002), and is well known in different contexts (Kroner 


1977). 


To obtain a formula for Y notice that (2.6) and (2.8) imply that 

0 = F 2 E^ = F 2 E 1 + F 2 E 2 = F 2 E 1 + F 2 L ^Fl 2 F 2 J^ (4-9) 

where the inverse of L is to be taken on the subspace 77. Solving for J' gives 


J' = -(F 2 L-in 2 F 2 )-'F 2 Ei, 


(4.10) 


where the inverse is to be taken on the subspace J. Then by applying IIi to both sides 
of this equation and equating IIiJ' = Ji with — YEi we obtain the desired formula 


Y = niF 2 (F 2 L-'n 2 F 2 )-'F 2 ni, 


(4.11) 


for Y, as given in formula (19.29) of Milton (2002[ ). 

Another formula for Y is obtained by taking an arbitrary constant zo 7 ^ 0, and 
dehning 

P' = J' - zoE'. (4.12) 


Applying Fi to both sides of (4.12) gives 


FiP' — —ZqEI — —2;o(Ei + E 2 ), 


( 4 , 13 ) 
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and applying 112 to both sides of (4.13) gives 


112?^ — J 2 — Z 0 E 2 ~ (L — zoI)E2. (4.14) 

Combining these results we see that P' satishes 

[Pi + zo{L - zoI)^'ni]P' = -zoEi. (4.15) 

Assuming that the operator [Pi + zo(L — 2;oI)^^PIi] is nonsingular this gives 

P' = -zo[Ti + zo{L - zoI)^'ni]-'Ei. (4.16) 

Applying PIi = I — PI 2 to both sides yields 

Ji - zoEi = -(Y + zoI)Ei = -zoTi[T, + zo{L - zol)-^U,]-^E, (4.17) 

As this holds for all Ei G V we obtain the formula 

Y = —zqPIi + ^ori[Pi + 2 :o(L — ^ol) ^IIi] ^PIi (4-18) 


which is a special case of the formula (19.37) obtained in Section 19.5 of Milton (2002). 


5 Multiplying superfunctions 

Multiplying superfunctions is similar the way electrical circuits, each with 2m terminal 
can be combined. An example is shown in Figure 

Suppose we have two superfunctions, (E®)' and (E^)": 

1C' = £'(SJ' = {V^y © (V°)' © H' with -H' = E( © E' © • • • © E', 

/C" = S" © J" = (V^" © (V®)" © E" with E" = E" © E^' © • • • © E(;, (5.1) 

where the spaces (V'^)', (V*^)', (V'^)", (V*^)" all have the same dimension m. To take their 
product one needs to hrst hnd two nonsingular linear operators and which map 
(V*^)' to (V'^)". The resulting product superfunction 

E^ = (E*)'xm(E^)", (5.2) 


is the subspace collection 


)C = g(Bj= (V^)' © (V°)" © E, 

(5.3) 

where 


E = E( © E2 © • • • © E) © E" © E2 © • • • © Ef, 

(5.4) 

and the operator L acting on E is 


j ^ 

(5.5) 


i=l £=l 
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Figure 5: Multiplying superfunctions is like hooking networks, with an equal number 
of input and output terminals, together in series. Shown in (a) and (b) are 6 terminal 
electrical networks, each (along with their respective tree-like battery conhgurations on 
the opposite side of the circuit board that are not shown here) represent a superfunction 
as the terminals have been divided into input terminals (/(, J^, and Jg for the circuit 
(a), and I 2 , and Jg for the circuit (b)) and output terminals ( 0 (, O 2 , and O 3 for 
the circuit (a), and O”, O 2 , and O 3 for the circuit (b)). The product superfunction is 
the 6 terminal electrical network (along with its tree-like battery conhgurations on the 
opposite side of the circuit board ) shown in (c). Note there is some hexibility in how one 
multiplies superfunctions: instead of connecting the terminals O' with J" for i = 1,2, 3, 
one could for example, connect 0 ^, 02 , and O 3 with any permutation of and Jg. 

This is why, when taking a product, one needs to specify the maps (M® and M"') one is 
using between the output space of one superfunction, and the input space of the second 
superfunction by which one is multiplying it. 


in which A' and A” are the projections onto P' and V'/. A vector E is in T if and only 
if we can hnd vectors 


E' = (E^y + (E°)' + E '2 G S', 

E" = (E^)" + (E^)" + E" e S", (5.6) 

such that 

(E^)" = M^(E^)', E = (E^)' + (E^)" + E '2 + E' 2 ', (5.7) 

with 


wy G (VO', (E")' G (V")', E '2 G n', (EO" G (VO", (E^)" G (V")", E' 2 ' G n". 

(5.8) 
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A vector J is in if and only if we can find vectors 


J' = (JO' + (J°)' + J'a e J', 

r = (J^" + (J^)" + J" G J", (5.9) 

snch that 

(J^)" = J = (J^)' + (J°)" + J '2 + J 2 , (5.10) 

with 


(jO' g (vO', (J"")' G (V^^)', J' G n\ (jO" g (vO", 


(jO)" e (yO)", 


J' 2 ' G W. 


(5.11) 

To ensure that the two spaces £ and J are independent we need to make the technical 
assumption that M® and are chosen so that the operator A mapping (V*^)' to (V'^)", 
dehned by 


A = M^(n^)'r; - (n^)"r;'[M^(n°)'r; + M^(n^)'r'2], (5.12) 


is nonsingular (i.e. the null-space of the operator contains only the zero vector). Our 
aim is to show that if A is nonsingular and 

E = (E^)' + (E°)" + E '2 + E' 2 ' = j = (J^)' + (J^)" + J' 2 +J' 2 ', with E G 3 ej (5.13) 


then E = J = 0. First note that by resolving (5.13) into components in the spaces (V^)', 
(V^)", Tf', and H” we obtain 


EO' = (JO', (E"")" = (J"")", E '2 = J' 2 , E'' = J 


r" 

’ 2 - 


(5.14) 


Also since E G and J G 57” there exist vectors (E*^)', (J*^)' G (V*^)' and (E^)", (J^)" G 

(V^)" such that (5.6) and (5.9) hold. Since £' H J' = {0} and £” fl J” = {0} it follows 

that 

P = (E^)' - (J°)' = E' - J' ^ 0 or E' = J' = 0, (5.15) 

and 

Q = (EO" - (JO" = E" - J" ^ 0 or E" = J" = 0. (5.16) 

Now we have 


(n^)'r;p = (n^)'E' = (e^)', (n^)'r' 2 P = -(n^)'j' = -(j 
(n0"r'i'Q = (n0"E" = (eO", 


rON/ 


(nO"r' 2 Q = -(n0"J" = -(J )"• (5.i7) 

Since (E^" = M®(E0)' and (J^" = M-^(jO)' we get from the hrst pair of equations in 

r;p, (jO" = -M^(n°)'r' 2 P, 


(5.17) the result that 

(EO" = 


which implies 


Q = [M^(n°)'r'i + M0n^)'r'2]P. 


1 


(5.18) 


(5.19) 


Substituting this back in the second pair of equations in (5.17), and using (5.18), gives 

(n0"r'i'[M^(n°)'r; + M^(n°)'r'2]p = M®(n°)'r'^p 


(nO"r'2[M^'(n^)'r'i + M0n°)'r'2]P = M^(n^)'r'2P. 


r05/-n/- 


(5.20) 
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These two equations are not independent since by adding them we obtain 


[M®(n°)T; + M‘^(n^)T'2]p = M®(n^)T;p + M-^(n°)T'2P 


(5.21) 


which is obviously satished. Also the hrst equation in (5.20) says P is in the null space 


of A, which by our assumption implies P = 0. Then (5.19) implies Q = 0 and this rules 


out the hrst possibilities in (5.15) and (5.16), implying E' = J' = 0 and E" = J" = 0. 
We conclude that E = J = 0. 

To check that the space £ ® J spans (V^)' © H , we just need to count 

dimensions. The dimension of the space on the right is 2m+dim('H). The dimension of 
£ according to (5.6) is dim(£^')+dim(T") less m because of the m constraints (E^)" = 
M®(E‘^)'. Similarly the dimension of £[ is dim(j7')+dim(j7’^')-m. Adding these up, we 
get the dimension of £’©^7 is dim/C'+dim/C"-2m=2m+dim('H')+dim('H')=2m+dim('H). 

Let F' and F' be the functions associated with the superfunctions {F^)' and {F^)". 
Given operators 

= L" = ^j"A'', (5 .22) 

i=l i=l 

where A' projects onto V[ and A'/ projects onto P", and given input helds (E^)' and 
(J^)' we can calculate 


(J^yj 1,(J7 

(E^)' = M^(E^)", 

_ T,. aeo" 
■J®)" J I 


(J^' = M-^(J 


Jl 


(5.23) 


From the knowledge of (E*^)' and (E^)', and of (E*^)" and (E^)", we can calculate the 


helds E', E", J', and J" of the form (5.6) and (5.9) solving the Y' problem and the W 
problem; 


E' G £', J' G j\ j; = l'e;, 

E" G £Y J" G J\ J" = L"E". 


(5.24) 


Then the helds E and J given by (5.7) and (5.10) solve the Y problem in the space /C, 


and the function associated to the superfunction F® is given by the product rule 

'M® 0 


F = F' 


0 M 


(5.25) 


Let us choose a basis (v{)", (v^)",... (v^)" for (V^)", choose a basis (vf)", (v^)",... (v(^)" 
for (V'^)", take M-®(vy)", M®(vy)",... M.^{-v^Y as our basis for (V’^)', and choose a ba¬ 


sis V 


oy 

1 ) ) 


rOy 


rOy 


for (V*^)'. Then the operator is represented as the identity 

matrix in the basis. Let us also choose the operator so it is represented by minus 


the identity matrix in this basis. Then in this basis the relation (5.25) takes the form 


F = F' 


I 0 

0 -I 


(5.26) 
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Note that we could have avoided this slightly awkward multiplication rule if we had 
replaced the dehnition (2.17) of the associated function by 


EO 

-jO 


= F 




(5.27) 


Then the multiplication rule (with this choice of M® and M"^) would have become simply 
F = F'F". We chose not to do this in the interest of preserving the symmetric roles of 
the spaces £ and J in the dehnition of the function associated with the superfunction. 

In passing, let us suppose there is an inner product on the vector spaces 1C' and /C”, and 
that the sets of spaces {S', J'}, {(V')', (V°)', ... ,V{}, [S", J"}, {(V')'', (V°)", V'{, V'i, 

all contain mutually orthogonal spaces. For any two helds 


p = + pO + p' + p", Q = + qO + Q/ + Q" 

in the vector space /C, with 

G (VO', p^, e P', Q' G W, P", Q" g U" 

let us dehne the inner product of them to be 

(P, Q) = (P^ Q^y + (P°, Q®)" + (P', Q')' + (P", Q")", 


(5.28) 


(5.29) 


(5.30) 


((J®)',E')'- ((J0",E")" 


in which ( , )' and ( , )" denote the inner product on the spaces 1C' and 1C" respectively. It 
is immediately clear from this dehnition that the subspaces (V^)', (V'^)", "PJ, V' 2 ^ ■ ■ .,P), 
V'{, V 2 , ■ ■ ■, Vy are mutually orthogonal in the new superfunction. Now take a held 
E G T and J G JT”. By the dehnition of these subspaces there must exist helds E' G £' 
and E" G £" such that (5.6) to (5.8) hold, and helds J' G J', J" G J" such that (5.9) to 
(5.11) hold. Consequently we have 

(J, E) = (J' + J" - (J®)' - (J^)", E' + E" - (E°)' - (E^)") 

= ((J^)', (E^)')' + ((J^", (E^")" - (J', (E°)')' - (J", (E")")" 

= -((J®)', (E®)')' - ((JO", (EO")" 

= -((J®)', (E®)')' - (M'^(J^)', M^(E°)')" 

= -((J®)', (E®)')' - ((M^)^M-^(J^)', (E^)')', 

in which (M®)1 is the adjoint of M®. So we see that the spaces J and £ will be 
orthogonal if we choose 

(M®)'^M^ = -I. (5.32) 

Note that the orthogonality of the spaces J and £ immediately implies that they have 
no nonzero vector in their intersection. 

In the case of nonorthogonal subspace collections, we are free to choose the maps 
M® and that map (V*^)' to (V'^)", so long as they and the map A are nonsingular. 
However, in view of (5.32), it would be quite natural to restrict our dehnition of multipli¬ 
cation by requiring that = —M®, i.e. one can pick a nonsingular map M mapping 
(V*^)' to (V^)" and set 


(5.31) 


= M, = -M. 


(5.33) 
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With this choice, subtracting the equations in (5.20) gives 


- r")M(n^)'(r; - r2)P = mp 


(5.34) 


Returning to the case where the subspaces are orthogonal, (5.32) is satished if MMl = 
—I. An alternative way to see that J and E have no nonzero vector in their intersection 
is as follows. Choose an orthonormal basis (vf)', (v^)',... (v^)' for and take 

as a map such that M®(vf)', M'®(v^)',... M®(v^)' form an orthonormal 
basis for (V^)". Then the operator M'® is represented as the identity matrix in the basis, 
andM*^ is represented by —I. Now, recalling the dehnition of the norm |Q| = (Q, 
of a vector Q recall that the action of the operators (H^)" cannot increase the 

norm of a vector, while P^^ — Fg and F'^^ — F 2 preserve the norm (as can be seen if we 
take a basis where these are diagonal). Hence (5.34) can be satished only when there is 
a P G such that 

(r; - f'2)p g (v°)' (F'l' - r'2')MP g (v 

Then as F) + F 2 = I and F'/ + F 2 = I we obtain 

(r; + f'2)p g (v°)' (f)' + f'2')mp g (v 


Adding and substracting (5.35) and (5.36) then implies 


F;P G (V")', F'2P G (V 


0\i 


ii\" 


{vr. 

( 5 . 35 ) 

{vr. 

( 5 . 36 ) 

F"P G (V^" 

( 5 . 37 ) 


which is excluded by our assumption that V has no vector in common with £' or J' and 
that V" has no vector in common with £" or J”. 


6 Multiplicative identity superfunctions 

Suppose we are given nonsingular maps and which map the m-dimensional space 
to the m-dimensional space Let /C" denote the 2m-dimensional space 

/C" = (V°)'© (VO". (6.1) 

Within this space dehne £" as the subspace consisting of all vectors of the form E = 
V + (M^)“^v with V G (V^)'' and dehne J" as the subspace consisting of all vectors of 
the form J = w + (M'^)“^w with w G (V^)". If these subspaces have a vector in common 
then 

V + (M^)“^v = w + (M'^)“^w, i.e., v — w = (M'^)“^w — (M'®)“^v. (6.2) 

In this last equation the helds on the left and on the right lie respectively in (V^)^' and 
(V*^)'. As the intersection of these subspaces consists of only the zero vector, we conclude 
that both sides must be zero, i.e. w = v and 

u == (M®)-^v = (M-^)-^v (6.3) 
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Thus, M®u = V = M-'u and if we assume that is nonsingular, then 0 = 

u = V = w. So under this assumption the subspaces have only the zero vector in their 
intersection. Then, since they each have dimension m we conclude that 

/C" = (V^)' © (VO" = £" © J", (6.4) 

which dehnes a superfunction (F^)" in which Ti is empty. 

We now look at the associated superfunction problem. As the space "H is empty, if 
we are given vectors and in the input space (V'^)" the superfunction problem then 
consists of hnding vectors and in the output space (V'^)' such that 

E^ + E® e T", e J”. (6.5) 


From our dehnition of the subspaces £” and J" we immediately see that the superfunction 
problem is solved with helds 


E“ = (M^')-^E', = (M-^)-^J^ 


implying, through (2.17), that the associated function is 

(M^)-i 0 


F" = 


0 M' 


r-l 


( 6 . 6 ) 


(6.7) 


So if we take another superfunction (F*)' and multiply it by this superfunction (F^)", 
the product rule (5.25) implies that the resulting superfunction has the associated 
function 

F = F'. (6.8) 


We conclude that this superfunction (F^)" is the multiplicative identity, when multipli¬ 
cation is dehned with the maps M'® and M'^. 


7 Addition of y-subspace collections and embeddings 

Adding superfunctions is similar the way electrical circuits, each with n terminals can be 
combined. An example is shown in Figure 

Suppose we have Y{j) and Y{k) subspace collections: 

/C' = £' © d7' = V' © H' with 77' = F( © F' © • • • © F', 

1C" = E" © J" = V" © F" with F" = F(' © F 2 ' © • • • © V'i (7.1) 

where the spaces V' and V" have the same dimension n. To dehne the sum of the 
subspace collections we need to introduce another n-dimensional space V and nonsingular 
operators S' and S" which respectively map V' and V" to V. Then the sum of the subspace 
collections 

/C = /C' +{s^s"} /C" (7.2) 

is the subspace collection 

/C = ^©j7 = V©F, (7.3) 
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CjZj 



(a) (b) 

^ 6^2 



(c) 


Figure 6: Adding F-subspace collections is like hooking networks together in parallel. 
The 4 terminal networks in (a) and (b), each representing (along with their respective 
tree-like battery conhgurations on the opposite side of the circuit board that are not 
shown here) Y (3) and Y (2) subspace collections, are added together to form the 4 ter¬ 
minal network in (c) which is a T(4) subspace collection. Note that the circuit in (b) 
is really only a 3 terminal network, so it has been embedded in a 4 terminal network 
(with no electrical connections to the 4th terminal). Also note there is some flexibility in 
how one adds together these subspace collections: we connected the terminals A', B', C, 
and D', to respectively the terminals A", B", C", and D'\ but we could have connected 
them to any permutation of these terminals. This flexibility is reflected in the need to 
introduce nonsingular operators S' and S" which respectively map V' and V" to V, before 
addition can dehned. 


where 

H © P' © • • • © P' © P" © P" © • • • © Vl (7.4) 

Here a held E = Ei + E 2 , with Ei G V and E 2 G B, is in S if and only if there exist 

helds 

E' = E; + E '2 G S', E" = E” + E” G S”, (7.5) 

with 

e; g V', E '2 G n', E)' G V", E' 2 ' G n", ( 7 . 6 ) 

such that 

s'e; = s"e;' = Ei. ( 7 . 7 ) 

Also a held J = Ji + J 2 , with Ji G V and J 2 G "H, is in 77 if and only if there exist helds 

J' = j; + J '2 G S', J" = J'l + J' 2 ' G S", (7.8) 
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with 

J'l 

G V', J '2 G J'l' G V", J' 2 ' G 

(7.9) 

such that 

S'J'i + S"J'i' = Ji. 

(7.10) 

So given Ei G V, we let Ej = (S') ^E and E'/ = (S") ^Ei, and we solve the Y-problem in 
each of the two subspace collections Y(j) and Y(fc), finding fields satisfying (7.5), (7.6), 
(7.8), and (7.9) with 

J '2 = L'E2, J '2 = L"E' 2 ', (7.11) 

where 




i=l i=l 

(7.12) 


and A' projects onto V[ while A'/ projects onto V”. Hence we have 


J 2 = J '2 +J" = L(E 2 + E"), with L = L' + L". 


Then (7.10) implies 


Ji = s'j; + s"j;' = s'y'e; + s"y"e" = yei, 


(7.13) 


(7.14) 


where 

Y = S'Y'(S')“^ + S"Y"(S")"^ (7.15) 

If we have a basis vi, V 2 ,..., v„ for V, then it is natnral to take (S')“^vi, (S')“^V 2 , ..., 
(S')~^v„ as a basis for V', and to take (S")“^vi, (S")“^V 2 , ..., (S")“^v„ as a basis for V". 
Then the operators S' and S" are represented by identity matrices, and in these bases 
(7.15) becomes Y = Y' + Y". 

In the case where either or both of the snbspaces V' and V" have dimension less than 
the dimension n of the subspace V we can first do an embedding. For example suppose 
V' has dimension n' < n. Then let W' be a space of dimension n — n'. Construct the 
subspace collection _ _ _ 

K,' = £' ®J' = V ®'H\ (7.16) 


where 


V' = V' © W', S' = £' ® W'. 


(7.17) 


Then given a field Ej G V we can express it as a sum Ej + W' with Ej G V' and 
W' G W'. We write E'j^ = TE'j^ where T is the projection onto V'. Given this Ej and 
solving the Y-problem associated with K,' we obtain helds E' and J' satisfying 


E' = e; + E'2 G £', E'l G V', E'2 G 

J' = j; + J'2 G J', J'l G V', J'2 = LE 2 G H. (7.18) 

It follows that the Y-problem in the space /C' is solved with fields 

E' = W + E' = W + E'l + E'2, and J' = J) + J'2 with J'2 = LE 2 , (7.19) 


implying that 


J'l = -ye; = -yte;. 


(7.20) 
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We conclude that the new F-problem has an operator Y = Y\&, i.e. its range is not 
the whole space V but only at most the subspace V'. After making such embeddings to 
ensure that V and V" (or rather V' and V have the same dimension as the dimension n 
of the subspace V, we are then free to add them. 

The additive zero is easy to find. Let us consider the degenerate subspace collection 

/C" = 8 " = V (7.21) 

Clearly 'H" contains only the zero vector, and we are forced to choose L" = 0. Given 
El G V". The Y-problem is solved with vectors 

E" = El, El = Ji = J 2 = J = 0. (7.22) 

Implying the associated Y-operator Y is zero: thus the subspace collection ( |7.21[ ) is the 
additive zero. Note that this subspace collection does not satisfy the property 8 " nv" = o 
which is needed for the inverse of Y to exist, which is not surprising since Y = 0 has no 
inverse. 

Now suppose we have a subspace collection 

]C' = 8 ' ®J' = V ®H' with W = V[®V'^®---®V-, (7.23) 

with associated operator Y'(2:(, • • •, when 

L' = ^^'A'. (7.24) 

i=l 

It is clear that if we replace L' by 

= (7.25) 

i=l 

then the solution to the Y-problem will give the Y-operator 

Y(-/„ -4,.... -4) = -Y(4,4,... ,4), (7.26) 


where to obtain this last identity we have used the homogeneity of the function. Since 
adding (7.26) to the associated operator z'^ we started with gives zero, it 

is tempting to conclude that we have found the additive inverse. However the function 
(7.26) is not the Y-operator valued function oi z'^^ z' 2 ^^ z'^ associated with the subspace 
collection (7.23), whose dehnition does not allow us to choose L' of the form (7.25). This 
is made more clear in the case where we have an orthogonal subspace collection since then 
the imaginary part of (V, Y(4(, ^ 2 , • • •, z'^N^ is generally positive when ^ 2 ) ■ ■ ■ iz'n ah 
have positive imaginary parts, and —Y(4(,Z2) • • • )^n) then does not share this Herglotz 
property. So the additive inverse of an orthogonal subspace collection should typically 
not be an orthogonal subspace collection. We will find the proper additive inverse in 
section [ 
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(b) -V\AA 

Zi 


kiZ4 —vw— 
gets replaced by ^W\A< ► < r 

-VW- 

^324 



Figure 7: Substitution of F- and Z-subspace collections is like replacing all resistors of 
one type by a compound network. If one takes a subspace collection, as, for example, 
represented by the 4-terminal network in (a) and replaces Zi by the network in (b), where 
ki + {l/k 2 + l//c 3 )“^ = 1, to ensure this replacement does effect the resistance when 
zi = Z 4 ^ = Z 5 = 1, one obtains the subspace collection as represented by the 4-terminal 
network in (c). 


8 Substitution of subspace collections 

Another familiar operation that we can do with rational functions is to make substitu¬ 
tions. Substitution of one subspace collection in another is similar to the way it can be 
done in electrical circuits. An example is shown in hgure Thus if Y (zi, Z 2 , ■ ■ ■, Zn) is 
a m X m matrix-valued homogeneous function of degree one and Z'[z'l, z^, ■ ■ ■ 1 Zp) is a 
scalar-valued function, say normalized with 

Z'(l,l,...,l) = l, (8.1) 


then 

Y'\z[, 4 ) • • • > 4 ’ ^2, • • • , ^n) = Y(Z(4, 4, . . . , 4), ^2, • • • , Zn) (8.2) 

will be another m x m matrix-valued homogeneous function of degree one. What is the 
analogous operation on subspace collections? It is natural to expect there should be one, 
just as in a network of n types of resistors one can replace each resistor of type 1 with a 
network of p other resistors. 

Extending the treatment given in Section 29.1 of Milton (2002) let us suppose that 
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we are given a V (n)-subspace collection 


/C = ^ © = V © Pi © P 2 © • • • © (8.3) 

and a (3,p)-subspace collection 

P' = W' © ^' © = P( © P' © • • • © p;, (8.4) 

in which V is m-dimensional and W' is one-dimensional. Let Y(zi, Z 2 , ■ ■ ■, Zn) and Z'{z[, Z 2 , ■ ■ ■, Zp) 
denote the functions associated with these subspace collections. We take as our new 
(2, n + p)-subspace collection, 

/C" = S" © J" = V" © P" © P" © • • • © P", (8.5) 

where 

S" = J” = {J (8.6) 

and 

V" = V © U\ 

V" = Pi © P' for 1 < f < p, 

= Vi+i-p®U' for p + l<i<n + p—1. (8.7) 

in which © denotes the operation of taking the tensor product of two subspaces. Vectors 
in the space 

JC" = 8" © J" = (/C © W) © (Pi © {S' © J')) (8.8) 

spanned by these subspaces are represented as a pair [P, u'] added to a linear combination 
of pairs of the form [Pi, P'], where P G /C, u' G U\ Pi G Pi, and P' G © J'. 

Now define 

P = Pi © P2 © • • • © Pn, 

= P" © P" © • • • © P", (8.9) 

and suppose that we are given solutions to the equations 

n 

J 2 = ^ ^ ^jAjE 2 with El + E 2 G Ji + J 2 G 77, Ei, Ji G V, E 2 , J 2 G P, 

7=1 

n 

j' + J' = A'(e' + E') with e', j' G U\ E' G 8 \ J' G J\ 
j=i 

(8.10) 

where 

Zi = Z{z[,z' 2 ,... ,z'p), (8.11) 

while Aj and A'- are the projections onto Vi and P). Let us introduce 

P, = A,E 2 , P' =A'(e' + E'), (8.12) 
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and set 


z'l = z[ for 1 < -i < p, 




= Zi+i-p for 

p + 1 < i < n + p — 1. 

(8.13) 

Then, in the 

new subspace collection, the 

! vectors 



E'l' 

= [El, e'] e V", 

E" = [E 2 , e'] + [Pi, E'], 



J'l' 

= [Ji, e'] G V", 

J" = [J 2 , e'] + [Pi, J'] 

(8.14) 

satisfy 







E" + E" G S 

", J'l' + J' 2 ' G J". 

(8.15) 

Additionally, 

we have 





n 


p n+p—1 



E2 = lE 

P., e'] - [Pi, e'] 

+ 1P,, Epa= E 

(8.16) 


2 = 1 2=1 2=1 


where 


P"= [Pi, P'] for l<i<p, 

= [Pj+i-p, e'] for p + l<i<n + p- l 


(8.17) 


satisfies P'' G V”. Similarly, and using the fact implied by (8.11) that j' = Ze' = zie\ 
we have 


n p n+p— 1 

3 « = 1^ ZjPi, e'] - [P,. j'] + [P,. j] 4P'] = zfp" e n". (8.18) 


2=1 


2=1 


2=1 


Given a basis vi, V 2 ,..., for V and a vector u' in W' it is natural to take (vi, u'),(v 2 , u'), 
..., (vm, u') as our basis for V". Choosing e' so that e' = u', it is evident that 
Y {Z’{z'i, Z 2 ,..., zp, Z 2 ,. ■ ■, Zn) is the matrix-valued function associated the new subspace 
collection, represented in these bases. 

There is a similar subspace operation corresponding to substituting the Z-function 
Z'{z\^ Z2,... ,Zp) into another Z-function Z(zi, Z2,..., Zn) to obtain 

Z"(z(, 4, . . . , Zp, Z 2 , . . . , Zn) = Z(Z(z[, Z 2 ,...,Zp),Z 2 ,..., Zn). (8.19) 

Given a Z(n)-subspace collection 


l~i — lA ® S ® J' — T’l © 1^2 © • • • © 'Pni 


( 8 . 20 ) 


and a (3,p)-subspace collection 

n' = U' ®J' = V[®V' 2 ®---®V'p, (8.21) 

in which U is m-dimensional and U' is one-dimensional, we take as our new (3, n + p — 1)- 
subspace collection. 


/C" = W" © S" © J" = © • • • © P", 


( 8 . 22 ) 



where 


U'' = U®U\ E” = {E®U')®{Vx®E'), J” = {J ®U')®{Vx®J'), (8.23) 

and 

V'l = Vi® V[ for 1 < i < p, 

= Vi+i-p®U' for p + l<i<n + p—1. (8.24) 

Suppose that we are given solutions to the equations 

n 

j -|- J = ^ ^ .ZjA.j(e -|- E) with e, j G W, E G E^ J G EJ•, 

i=\ 

n 

j' + J' = A'(e' + E') with e', j' G W, E' G E', J' G J', 

i=i 

(8.25) 

where Zi = ZZ 2 , ■ ■ ■ ^ , while A* and A'- are the projections onto Vi and V'y Let us 

introduce 

P, = A,(e + E), P' = A'(e' + E'), 

and dehne z” by (8.13), and P" G V” by (8.17). In the new subspace collection, the 
vectors 


satisfy 


e" = [e, e'] G W", E" = [E, e'] + [Pi, E'] G E" 
j" = [j, e']GW", J" = [J, e'] + [Pi, J']g:7" 


+ E" = [5^ Pi, e'] + |Pi, ^P'] - [Pi, e'l 


Z=1 

n 


3=^ 

P 


EPi, e'] + iPi. Eli 


i=2 

n+p—l 




E L. 


2=1 


and, using (8.11), 


(8.26) 


(8.27) 


j" + J" = lE e'] + |Pi. E ‘2LI - |Pi. j'l 


i=l 


i=i 

p 


E^iPi. e'] + [Pi, E^p;i 


i=2 

n+p—l 


J=1 


E LP?. 


2=1 


(8.28) 


Given a basis Ui, U 2 ,..., for 77 and a vector u' in 77' it is natural to take (ui, u'),(u 2 , u'), 
..., (um, u') as our basis for 77". Choosing e' so that e' = u', it is evident from (8.26) that 
Z(Z'(4,4,...,zL) , ^ 2 , • • •, Zn) is the matrix-valued function associated the new subspace 
collection, represented in these bases. 
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9 Some other elementary operations on subspace col¬ 
lections 

A further operation we can do on functions Y( 2 ;i, Z 2 , ■ ■ ■, Zn) while retaining the homo¬ 
geneity of degree 1 in the variables zi,Z 2 , ■■■ ,Zn is to replace the function by [Y(l/2;i, 1 / 2 ^ 2 , ■ 
The analogous operation on the associated Y (n)-subspace collection is to interchange the 
subspaces S and J'. Similarly in a Z(n) subspace collection, interchanging the subspaces 
£ and JT” corresponds to replacing Z(zi, Z 2 , ■ ■ ■, Zn) by [Z(l/zi, 1/^2, •••, as noted 


in Section 29.1 of Milton (2002). We call such a transformation a duality transformation. 

ula 

( 9 . 1 ) 


As a consequence of the duality transformation (4.3) immediately implies the formula 

ri)|-'r„. 


z-' = r„[(r„ + r,)L(rc 

One simple thing we can do in a function is set Zj = Zk'- the analogous operation in a 
subspace collection is to replace Vj Q)Vk hy a. single subspace. 

Another elementary operation we can do on a Z{n) subspace collection is as follows. 
Let U be expressed as the direct sum 


U = U'®W, 

which dehnes the projection $ onto W. We now take as our subspace collection 

n = u' ®£®j' = Vi®V2®---®Vn, 

where 

y' = y © W. 


( 9 . 2 ) 


( 9 . 3 ) 


( 9 . 4 ) 


Then any solution to the Z-problem (2.23) with e immediately generates a solution 


to the Z-problem associated with the subspace collection (9.3): 



j'eW', Ee^, J'ey', j'+ J' = L(e + E), 

( 9 . 5 ) 

where 

n 



L = j' = ^j, J' = J + (I-$)j, 

( 9 . 6 ) 

which ensures that 

j + J = j' + J' and (I - $)j G W. 

Hence the new subspace collection has a Z-operator 

( 9 . 7 ) 


Z' = $Z, 

( 9 . 8 ) 

when applied to helds in W. 
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Realizing any y-matrix with elements 
rational functions of degree 1 

that are 

Given any homogeneous rational function of degree 1, 



4 p{Zi,Z2,...,Zn) 

Z{Zi,Z2,...,Zn) = , 

q{Zi,Z2, ...,Zn) 

(10.1) 


, 1/^n)] 


-1 
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satisfying the normalization Z’(l, 1,..., 1) = 1 where p(^i, 2 : 2 , • • •, Zn) and q{zi, Z 2 , ■ ■ ■, Zn) 
are homogeneous polynomials of degree k and k — 1 respectively, where fc is a positive 
integer, our hrst goal is to hnd a Z{n) subspace collection 

H = hi ® S ® = 'Pi © P 2 © ■ ■ ■ © Pm (10.2) 

where U is one-dimensional which has Z{zi, ^ 2 , • • •, Zn) as its associated function. Without 
loss of generality we could set Zn =1, and thenp( 2 ;i, 2 : 2 , ■ ■ ■, Zn-i, 1) and q{zi, 2 ^ 2 ,..., Zn-i^ 1) 
are just arbitrary polynomials of the n — 1 variables 2 ^ 1 , 2 : 2 ,..., Zn-i- Also without loss of 
generality we can assume 


(3.25) implies we can realize 


p(l,l,...,l) = g(l,l,...,l) = l. (10.3) 

The hrst step is to realize Z( 2 :i, 2 : 2 ,1) = 2 : 1^2 as an associated Z-function. Note that 

Z(zi,l) = zl (10.4) 

(10.5) 


and (3.26) implies we can realize 

^(2^1, 2:2) = C 2 ;i + (1 - c)^ 2 , 

for any constant c. Hence, by substitution we can realize 

Z(zi, Z2, 1) = 9(22)1/3 + Z 2 / 3)^/8 — (2^1 — ^2)^/8 = ^1^2- 
Making further substitutions, we can realize any product of the variables 

Z(2;i, 2)2, ... , Zn-l, 1) = ^i'2:2" . . . 


( 10 . 6 ) 


(10.7) 


where the a* are nonnegative integers. By repeated substitution in (10.5) we can realize 
any linear combination of such terms with coefficients summing to 1. Thus we can realize 
the polynomials p{zi,Z 2 ,..., ^n-i, 1) and g( 2 :i, 2 : 2 ,..., 2 :„_i, 1). 

Furthermore (3.25), with the roles of 2:1 and Z 2 interchanged, implies we can realize 

Z(^i,l) = 1/^1, (10.8) 


which by substitution into (10.6) implies we can realize 

Z(^i,2;2, 1) = Z 2 IZ 1 . 


(10.9) 


Substituting p{zi, Z 2 , ■ ■ ■, Zn-i, 1) for 2:2 and q{zi, Z 2 , ■ ■ m Zn-i, 1) for 2:1 we see we can hnd 
a subspace collection which realizes 


Z(2:i, 2:2, . . . ,2:„_i, 1) = 


p(2;i,2;2,...,2;„_i,1) 


( 10 . 10 ) 


g(2;i,2:2,...,^n-i,l) 

as its associated Z-function when Zn = 1. When Zn is not 1, the subspace collection will 
by homogeneity realize the function (10.1). 


Now from (3.13) we can realize 


Y(.,) = «) . 


( 10 . 11 ) 
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and realize 


0 ai2^2 
0 0 


( 10 . 12 ) 


Y(;^2) = 

By substitution of subspace collections, we can realize any ^-matrix where in the above 
formulae Zi and Z 2 are replaced by any normalized rational homogeneous functions of 
degree 1 (normalized in the sense that they take the value 1 when all variables take the 
value 1). Finally, by making suitable embeddings and adding subspace collections we can 
realize any F-matrix with elements that are homogeneous rational functions of degree 1: 
(10.11) with the appropriate substitutions realizes each diagonal element, while (10.12) 
with the appropriate substitutions realizes each off-diagonal element. 


11 Extension operations on subspace collections 

Let us suppose we have a Z(n) subspace collection 

'hi = lA (B £ ® = 'Pi © P2 © • • • © Pm 

where lA is m-dimensional. Let V be another m-dimensional space, and consider the 
space 

/C = V©H. (11.2) 

Suppose there is a nonsingular mapping T from lA to V. Define the subspace £ to consist 
of all vectors spanned by u + Tu as u varies in lA. Define J to consist of all vectors 
spanned by u — Tu as u varies in lA. Clearly we have 

V@U = £@J, (11.3) 

and consequently we obtain the Y (n) subspace collection 

P = £' ®J' = V®Pi®p2®---®Pm ( 11 - 4 ) 

in which 

£' = £®£, J' = J®J. (11.5) 

Furthermore given vectors satisfying 

j + J = L(e + E), E G J G JL, e, j G W, (11.6) 

where 

n 

j = Ze, 'L = ^ZiXi, (11.7) 

£=i 

we can set 

E 2 = e + EG?7, Ei = Te GV, J 2 =j + JG?7, Ji = —Tj. (11.8) 

Then we have 

El+ E 2 = Te + e + E G T', Ji + J 2 =-Tj + j + J G J7', (11.9) 
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and 


Ji = -YEi, with Y = TZT'h (11.10) 

Given a basis ui,U 2 ,..., for U, with respect to which the matrix valued function 
Z(^i, ^ 2 , • • •, Zn) is dehned, it is natural to take Tui, Tu 2 ,..., Tu^ as our basis for V. 
Then T is represented by the identity matrix, and the functions 7i{zi^ Z 2 , ■ ■ ■ ^ Zn) and 


Y(^ 1 , Z 2 ., ■ ■ ■ i Zn) are identical. We call the subspace collection (11.4) the extension of the 


subspace collection (11.1). 


12 Reference transformations and additive inverses 

Given the impedance network illustrated in Figure we are free the change the scaling 
constants Cj assigned to each bond to new constants c' and accordingly replace Zi with 
z[ = ZiCijc^ without changing the overall electrical response of the network. Analogously, 
given a homogeneous rational function Y {zi, Z 2 , ■ ■ ■, Zn) of degree one, an operation which 
preserves the homogeneity is obviously to multiply the variables by constants to obtain 
the function 

Y\z[, z^,..., z'J =Y{diz[,d2z'^,... ,dnz'J. (12.1) 

The associated operation on the Y (n) subspace collection (T, J') and (V, Vi, V 2 , • • •, Vn) 


is found by generalizing the analysis given after (29.3) in Milton (2002) and is as follows. 
Given nonzero (possibly complex) constants cf and c/, i = l,...,n we introduce the 
linear transformations 




niP + ^cfA.P, ^.■'(P) = n,P + 5^c/AiP, (12.2) 

i=l i=l 

on helds P G /C, where Ai is the projection onto Vi. These transformations leave the 
subspaces V and Vi invariant. Dehne the spaces 

£' = ^P^{£) and j' = ^\j). (12.3) 

These will have the same dimension as £ and J respectively. To see this, suppose 


'0'®(E) = for some E, E' G £. Then '0®(E — E') = 0 and since (12.2) implies 

= 0 only when P = 0 we conclude that E = E'. We need to make the technical 
assumption that 

^ V''^(J), for all nonzero E G J G J, (12.4) 

to ensure £' and J' have no nonzero vector in common. A more insightful meaning to 


the condition (12.4) is given in the next section. 


Let {£', J') and (V, Pi, V 2 , • • •, Vn) be our new subspace collection. Given a solution 
to the equations 


E G T, J G JT” , (I ~ Til) J — ^ ^ ^iAjE, 


(12.5) 


i=l 


in the original subspace collection, in which IIi is the projection onto V, the fields 
E' = ?/^®(E) and J' = '^‘^(J) will be a solution to the equations 


E'ef', J'eJ', (I-ni)J' = ^j'AjE', 


( 12 . 6 ) 


2=1 
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in the new subspace collection with 


' J / 

= ZiCi jc, 


(12.7) 


Since IIiE' = IIiE and IIiJ' = IIiJ, it follows that the Y-tensor functions of the two 


subspace collections are related by (12.1) where 

= cf/c/. 


( 12 . 8 ) 


In particular, if we choose cf = —c/ for all i we obtain di = —1. Then using the 
homogeneity of the function Y(zi, Z 2 , ■ ■ ■, Zn) we see that 


Y'(z;, 4,..., z'J = Y(-z[, -z', -4) = -Y(4,4,..., 4). 


(12.9) 


So if to another subspace collection, with an associated function Y"(zi,Z 2 , ■ ■ ■ ,Zn), we 
add this new subspace collection according to the prescription given in Section then it 
produces a subspace collection with an associated Y-function which is obtained by sub¬ 
tracting Y (zi, Z 2 ,..., Zn) from Y"(zi, Z 2 , ■ ■ ■, Zn)- In other words, when cf = —cj for all z, 
the subspace collection with subspaces (T', J') and (V, Vi,V 2 , ■ ■ ■, Vn) is the additive in¬ 
verse of the original subspace collection, having subspaces (T, J) and (V, Pi, V 2 , • • •, Vn), 


where and are linked by (12.3). If the technical condition (12.4) is not 


satisfied it appears that the subspace collection has no arithmetic inverse. 


13 


Operations on subspace collections leaving the as¬ 
sociated function invariant 


Note from (12.8) that if we choose c/ = cf for all i then the associated function remains 


invariant. More generally, if we are interested in leaving the associated function invariant, 
we could take 


£' = c£, j' = cj, v' = cv, n' = cn, v[ = cVi, (is.i) 

where C is a nonsingular linear operator which maps /C to itself. Then the helds E' = CE 
and J' = CJ will be a solution to the equations 

n 

E'e£', (I - n))J'= ^ ^,A'E', ( 13 . 2 ) 

i=l 

where 

n; = CHiC-f A' = CAiC-^ ( 13 . 3 ) 

are the projections onto V and V[. If vi, V 2 ,... is a basis for V then setting v' = Cvj 
we can take vj,V 2 ,.. as a basis for V. Since multiplying by C is a linear operation 
the coefficients in the expansions 

m m m m 

n;E' = 44 , HiE = n;J' = 44 , n.j = (i 3 . 4 ) 

2=1 2=1 2=1 2=1 
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can be equated: 


(13.6) 


and as a consequence the same matrix Y whose coefficients govern the relation 


Ji ^ ^ 


2=1 


also govern the relation 




(13.6) 


(13.7) 


2=1 


Due to this equivalence it suffices in the preceeding section to limit attention to the 
transformations ( 12 . 2 ) having c/ = 1 for all i: it is only the ratio di = cf jcf that has 
any real signihcance. Then '0'^(P) = P, and the technical condition (12.4) is violated 
only when there are nonzero vectors E G and J G such that 


E — El + E 2 , J + Ji + J 2 , El — Ji G V, J 2 — LE 2 G H, with L — ^ ^ 


2=1 




(13.8) 


Thus either Y(cf, cf,..., c®) has an eigenvalue of —1, or the E-problem with Zi = cf 
for all i has a nonunique solution (with a nontrivial solution having E 2 7 ^ 0 for the 
homogeneous problem with Ei = 0 and also Ji = 0 , the latter not being needed for 


nonuniqueness but being needed if (13.8) holds). If we are looking for the arithmetic 
inverse we take cf = —1 for all i, and the inverse exists except when Y(— 1 , — 1 ,..., — 1 ) = 
—Y(l, has eigenvalue —1 or when the Y-problem with Zi = 1 for all i has a 

nonunique solution (with the homogeneous problem having a nontrivial solution with 
both El and Ji being zero). 

There is a similar invariance of matrix functions associated with Z{n) subspace col¬ 
lections under the linear transformations, 


W = CU, 8' = C8, J' = CJ, V' = CVi. 


(13.9) 


These invariances are quite natural, as they are isomorphic to changing the basis in 
the vector spaces "H or /C. Thus, up to these trivial equivalences, the arithmetic inverse 
dehned in the previous section is unique. 


14 Multiplicative Inverses of superfunctions 

To hnd the multiplicative inverse of a superfunction {F^)' we let 1C" be a vector space 
with the same dimension as /C', and we take C as a nonsingular map from K,' to 1C". We 
then let 


J" = c{j'), n" = cn\ 

(V^)" = C(V^)', (V°)" = C(V0', Pf = CP'for i = 1,2,...J. (14.1) 

Introduce the transformation 

^(p) = n;'p - n'', (14.2) 
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where 11'/ is the projection onto (V'^)"© and 11/ is the projection onto 1-i" . This is a 
special case of the transformations in (12.2). We let £" = Note that the outpnt 

space (V*^)' gets mapped to the input space (V^)", and the input space (V^)' gets mapped 
to the output space and apart from these switchings we have essentially made an 

additive inverse in the F-problem. We still require the technical condition mentioned in 
the last section, to ensure that this additive inverse exists: the operator Y(l, 1,..., 1) — I 
is nonsingular and the Y-problem with Zi = 1 for all i has a unique solution (or more 
precisely the homogeneous problem does not have a nontrivial solution with both Ei and 
Ji being zero). 

Now suppose we are given a solution to the superfunction problem associated with 
(F*)', 

E' = (E')' + (E°)' + E '2 e J' = (J')' + (J°)' + j;, e J', J' = L'E' 2 , 


where 


L' = 4K 


(14.3) 

(14.4) 


2=1 


(14.5) 


in which A' is the projection onto P', and 

(E^', (jO' g (V')', (E^y, (j^y e (v^y, e/, j/ e n'. 

Now take vectors 

E" = y^iCE'), J" = -CJ', E'' = -CE/, J/ = -CJ/ 

(E^)" = C(E°)', (E°)" = C(E0', (J^)" = -C(J°)', (J^)" = -C(J^)'.(14.6) 

These solve the superfunction problem associated with (E®)", 


E" = (E^)" + (E°)" + E/ G £", J" = (JY + + J/ G J" 


fO\// 


where 




z" = 4 , 


J'' = L"E''. (14.7) 
(14.8) 


2 = 1 
1 // 


in which A- is the projection onto "P", and 


E')", (E)" G (V 


,i\// 


(E 


0 \// 


tO\// 


G(V 


0 \// 


E/, J/ G n". 


(14.9) 


Next let Ml denote the restriction of C to the subspace (V*^)', i.e., that op erator mapping 
(V*^)' to (V'^)", such that MiP = CP for all P G (V*^)'. Then from (14.6) we have 
(E^)" = Mi(E‘^)' and (J'^)" = —Mi(J‘^)'. To see that (E*)" is the inverse of the 
superfunction (E^)' when 




(14.10) 


2=1 


2=1 


we introduce the operator M 2 which is the restriction of to the subspace (V*^)", i.e., 
that operator mapping (V*^)" to (V^)', such that M 2 P = CP for all P G (V^)'. Then 
upon taking the product of the superfunctions (14.6) implies 

(Eoy'\ _ /{eY 
{joy I (j^)' 


(14.11) 
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where 


F = 


AM 2 )-' 0 \ 

V 0 -Ma’V 


(14.12) 


is the multiplicative identity operator. From this analysis it looks like there are many 
multiplicative inverses, paramerized by C, but in fact all are equivalent: this follows from 
the previous section. 


15 Pruning the subspace collections 


If an m terminal resistor network has a cluster of resistors which is not connected to the 
rest of the network, and that cluster does not have any terminals, only internal nodes, 
then clearly we can discard it without affecting the helds in the rest of the network and 
its response matrix. The analogous operation on subspace collections is called pruning. 

When L is close to zqI we can expand the inverses in (4.5) and (4.7) to obtain the 
series expansions 


E = 




(15.1) 


i=i 


Z = zoTo + 5^ro(L - 0 oI)[ri(L - zol)/zoYT^. (15.2) 

i=i 

From these expansions it is evident that is only those helds in "H that arise from products 
of the operators Fi, Ai, A 2 , ..., A„ applied to helds in U have any role in determining 
E and the associated function Z( 2 :i, 2 : 2 , • • •, Zn) (also j and J): so we may as well prune 
away any other helds from the vector space "H. Thus we can redehne "H as the smallest 
subspace containing U that is closed under the action of Fi, Ai, A 2 , ..., A„ and redehne 


£ = T{H, J = T 2 H, Vj = KjH, j = l,2 ,...,n. (15.3) 

This imposes constraints on the dimensions of these subspaces, as noted in Section 29.2 
of Milton (2002) where the results are given in the case where U has dimension 1 and 
where the spaces are orthogonal. Let pj be the dimension of Vj, j = 1, 2,... ,?7 ,, and let 
m, qi and q 2 represent the dimensions of U, £ and J. The total dimension of the vector 
space T-i is therefore 

h = m +qi +q 2 = P 1 +P 2 + ...+ Pn- (15.4) 


Now the space 

[Ai(7/ © £y\ © [A2(7/ © £y\ © ... © [A„(7f © £)\ (15.5) 

certainly contains W, and is closed under Fi (because it contains £) and is closed under 
Aj for each j. It therefore must be "H and A.j{U®£) which has at most dimension m + gi 
must be Vj. Therefore for each j we have the inequality 


Pj <m + qi, 


(15.6) 


and by summing these over j we see that 


q 2 <{n- l)(m + qi). 


(16.7) 
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Similarly the subspace 

[Ai(W © J)] © [A2(W © JT)] © ... © [A„(W © J)] (15.8) 

can also be identified with Ti and we obtain the inequalities 


Pj <m + q 2 , qi < {n - l)(m + ^ 2 )- 


(16.9) 


In the particular case when n = 2 the constraints (15.7) and (15.9) imply that the 
dimensions of the subspaces S and 77 can differ by at most m. Also in the case n = 2 we 
have 

Pi = (m + ga - P2) + qi = {m + qi- pa) + ^2 > max{gi, ga}, (15.10) 


and similarly for pa. 

Likewise we can redefine K. as the smallest subspace containing V that is closed under 
the action of Fi, Ai, Aa, ..., A„ and redefine 


^ = ri/C, 77 = r 2 /C, Vj = Aj}C, j = l,2,...,n. (15.11) 

Let V be the dimension of V, Pj be the dimension of Vj, j = 1,2,... ,n, and let gi and 
ga represent the dimensions of S and 77. The total dimension of the vector space /C is 
therefore 

h = gi + ga = n+pi ©Pa + ...+Pn. (15.12) 


The space 


V © [Ai(T)] © [A2(T)] © ... © [An,(T)] 


(15.13) 


certainly contains V, and is closed under Fi (because it contains S) and is closed under 
Aj for each j. It therefore must be /C and Aj{S) which has at most dimension gi must 
be Vj. Thus for each j we have the inequality 


Pj < 7 i, 


(15.14) 


and summing these over j we obtain 

g 2 < "^ + ('^ - l)gi. 


(15.15) 


Similarly since 


/C = V © [Ai(77)] © [A2(77)] © ... © [An{J)], 


we obtain the inequalities 


(15.16) 


Pj < g2, gi < n + (n - l)g 2 . 


(15.17) 


When n = 2 the constraints (15.15) and (15.17) imply that the dimensions of the sub¬ 
spaces S and 77 can differ by at most v. Also in the case n = 2 we have 


Pi = (72 -P 2 ) + qi-v = (gi -p 2 ) + q 2 -v> maxjgi, ga} - v, (15.18) 

with a similar inequality for pa. 
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16 Expressions for the numerator and denominator 
in the rational function 


Assume that a Z{n) subspace collection, with m = 1 has been pruned. Let wi,W 2 ,... ,Wgj+i 
be a basis ioTU®£ with wi in W and W 2 ,W 3 ,... ,Wqj+i in £. In this basis (ro + ri)Aj(ro + 
Ti) is represented by a (gi + 1 ) x (gi + 1 ) matrix Aj, and since the A* sum up to the 
identity operator it follows that 

n 

5;A, = I, (16.1) 

^=1 


Also, because the subspace is pruned, Aj(W ©£^) can be identified with Vi which implies 
the matrix Aj must have at most rank pi. It is exactly p* if "Pj fl JL = 0. The formula 
(9.1) for the Z-function implies 


1/Z(2;i, 2:2, • • • ,^n) = ei ■ [^ZjAj] ^ei, (16.2) 

i=l 


where ei is the qi + 1 component unit vector [1,0, 0,... 0]^. Hence, following the argument 
given in Section 29.2 of Milton (2002| ), Z{zi, Z 2 , ■ ■ ■, Zn) can be expressed in the form ( | 10 . 1 [ ) 
with numerator 


p{zi, Z 2 ,..., Zn) = det[^ ZiAi] = ^ (16.3) 

i=l ai,a2,...,an 

of degree 1 + qi, in which the sum extends over all oi, 02,..., On with 

n 

Oj = 1 + gi, 0 < Oj < Pi for i = 1, 2,..., n. (16.4) 

i=l 

Typically one expects that the maximum power of Zi in this polynomial will be the rank 
of Aj. However, for example, note that for the matrices 

/O 0 0\ 

Mi= 111 , M2 = I-Mi, (16.5) 

VO 11/ 

the maximum power of Zi in 

det[2;iMi + Z 2 M 2 ] = det[(2;i - Z2)Mi + Z2l\ = Z2[z2 + 2^2(© - 2 ^ 2 )] (16.6) 

is 1 while Mi has rank 2. 

Next let wi,Wgj+ 2 v • • ,^h be a basis for W © JT” with wi in U and Wg^_|_ 2 , Wgj+ 3 ... ,w/j 
in J'. In this basis (Fq + r 2 )Ai(ro + r 2 ) is represented by a (g 2 + 1) x (^2 + 1) matrix 
Bj, and since the Aj sum up to the identity operator it follows that 

n 

^B. = I, (16.7) 

i=l 
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Also, because the subspace is pruned, Aj(W © J) can be identified with Vi which implies 
the matrix Bj must have rank at most Pi. It is exactly Pi ii Vi Pi £ = 0. The formula 


(4.3) for the Z-function implies 

n 

Z{zi,Z2,.. . ,2„) ^ 02 ■ [© B,/^,]~^e2, 


(16.8) 


2=1 


where 62 is the q 2 + ^ component unit vector [1, 0,0,... 0]^. The denominator of this 
expression, as a polynomial in the variables 1 /zi, is 


det[^Bi/2;i] = ^ I3b,b2...b„/ 

i 1 


^■ 1^2 ■ ■ ■ 


(16.9) 


in which the sum extends over all bi,b 2 ,... ,bn with 

n 

^ 6 j = l + g 2 , ^<bi<Pi for i = 1 , 2 ,... ,n. 


(16.10) 


2=1 


Consequently, for the denominator in the expression (10.1) for Z( 2 :i, 2 : 2 , • • •, Zn)., we can 
make the identification 


q{zi, Z2,. . . , Zn) = /^6ib2...fen2:r ^'4" ^”1 

bi,b2,...,b„ 


(16.11) 


which is a polynomial of degree h — (1 + ^ 2 ) = <?i- Furthermore the identities (16.1) and 


(16.7) imply the polynomial p and q satisfy the normalization (10.3), i.e. 


E 


do 




^ ^ Pblb2...bn 1 - 

bl,b2,...,bri 


(16.12) 


17 The correspondence between rational functions 
of one variable and Z(2) subspace collections with 

m = 1 

In the case m = 1 and n = 2 there are two cases to consider. When the dimension of 


h is even, h = 2d, then in order to satisfy the inequalities (15.6), (15.7) and (15.9) the 


subspaces £ and 77 must have dimension d and d — 1 or vice versa and the subspaces 
Vi and V 2 must have dimension d. Without loss of generality, by making a duality 
transformation if necessary, let us suppose £ has dimension d — 1. Given u G W let us 
take as our basis for 77 the vectors 


V 2 ;-l = (riA,)." 'u, V 2 j = (Airi).' 'A,u, j = l,2,...,d, 


(17.1) 


SO that 


Vi = u, V2j = AiIZ 2 j_i, j = l,2,...,d, V2j+1 = riIZ2j_i, j = 1, 2,..., d - 1. 


(17.2) 
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These fields are independent since if they were not we could prune the subspace collection. 
The vectors V 2 j+i, j = 1,2,... ,d — 1, which number d — 1, must form a basis for S and 
so it follows that 

d-l 

r iV 2 d = ^ 7iV2i+i. (17.3) 

1 = 1 

Also we have 


roVi = vi, roV 2 j = (5jVi, j = 1,2,... ,ci, roV2j+i=0, j = 1,2,... ,d-l. (17.4) 

The 2d — l constants 71 ,..., ■jd-i and (5i,..., <5^ characterize the geometry of the subspace 
collection. The field e + E must have the expansion 


e + E — niV2j_i, 

i=l 

and consequently, setting Z 2 = 1 

d d 

j + j = [I + (^1 ~ l)-^l] (^ + E) = ajV2i-l + {zi — 1) CliV2i- 

i=l i=l 

Since ri(j + J) = 0 we arrive at the equations 

d d—1 d—1 

0 = ^ aiV2i-i + {zi - 1) ^ aiV2i+i + {zi - 1) ^ ad7iV2i+i 
i=2 i=l 

d-l 

= ^[Oi+i + ai{zi - 1 ) + liadizi - I)]v2i+1. 

implying 


(17.5) 


(17.6) 


2=1 


2=1 


(17.7) 

(17.8) 


Oj+i + ai{zi - 1) + 'littdizi - 1) = 0, i = l,...,d-l. 

Choosing a normalization with = (1 — these equations are solved with 

d-l 

ai = {l- zd'-^ - Y, Id-i+i-jil - ziY. (17.9) 


j=i 


Since 


we obtain 


2=1 


ro(e + E) — OiVi, ro(j + J) — [oi + {zi — 1) 

{zi - 1) Eti 


^(^ii 1) — 1 + 


Oi 


(17.10) 


(17.11) 


Conversely suppose we are given a rational function Z{zi, 1) with a denominator of degree 
at most d—1 and a numerator of degree at most d satisfying Z{1, 1) = 1. It can be 
expressed in the form 




ETMi-z.y^^ 


(17.12) 
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Comparing this with (17.11) we can make the identihcations 


d-l d-l 

1 - ZiY = ai = 1 - ^7d_j(l - zi)\ 

i=i i=i 

d—\ d 

= {zi - l)^(5iai 

j=0 i=l 

d—1 d—1 j 

^i'yd-l-^i—j {l-z,y+\{17.13) 

j=0 j=0 i=l 

which imply 


j 

Sj tg tj *^i+l ^ ^ ^i'^d—l+i—j j 1,..., h 1. (17.14) 

i=l 

Given the coefficients s and t we can inductively uniquely determine the coefficients 7 
and 6 : 


j 

7i ^d—jy ^0) ^j+1 ^ ^ i j 1,..., d 1. (17.15) 

i=l 


On the other hand when the dimension of h is odd, h = 2d — l, then in order to satisfy 
the inequalities (15.6), (15.7) and (15.9) the subspaces £ and J must have dimension 
d — 1 and the subspaces Vi and P2 must have dimension d — 1 and d or vice versa. 
Without loss of generality let us suppose Vi has dimension d — 1. Given u G W let us 
take as our basis for Ti the vectors 


V2J-1 = (r 1 Ai)^ ^u, j = 1, 2 ,..., d - 1, V2j 
which satisfy 


(AiFi)^ ^Aiu, j = 1,2 ,... ,d, 

(17.16) 


vi = u, V2j = AiJ/2j-i, V2j+i = riJ/2i-i, j = l,2,...,d-1. (17.17) 

Again these helds are independent since if they were not we could prune the subspace 
collection. The vectors V 2 j,j = 1, 2,..., d — 1, which number d—1, must form a basis 
for Vi and so it follows that 

d-l 

AiV2d_i = ^7iV2i. (17.18) 

i=l 

Also we have 


roVi = vi, roV2 j=djVi, j = 1 , 2 , ...,d- 1 , roV 2 j+i = 0 , j = l, 2 ,...,d- 1. 

(17.19) 

The 2d — 2 constants 71 ,..., ■jd-i and di,..., d^-i characterize the geometry of the sub¬ 
space collection. The held e -|- E has the expansion (17.5) and so, with Z 2 = 1, 

d 


d-l 


d-l 


j+J — [I+(^i~l)Ai](e-|-E) — aiV2i_i-|-(^i —1) aiV2i-|-(^i —1) arf7iV2i. (17.20) 


i=l 


2 = 1 


2=1 
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Since FiQ + J) = 0 we arrive at the equations 


d d—1 d—\ 

0 = ^ aiV2i-i + {zi - 1) X] aiV2i+i + {zi - 1) ^ ad7iV2i+i 
i=2 i=l i=l 

d-1 

= + ai{zi - 1) + 'Jittdizi - l)]v2i+i, 

i=l 


implying (17.8) which has the solution (17.9). Since 


(17.21) 


d-1 


d-1 


ro(e + E) = aivi, ro(j + J) = [ai + {zi - 1) + ad7i)]vi = [oi - 


2=1 


2=1 


we obtain 


Z{zi,l) = 1 - 


E?=i 

(l\ 


(17.22) 


(17.23) 


Conversely suppose we are given a rational function Z{zi, 1) with a denominator of degree 
at most d—1 and a numerator of degree at most d—1. It can be expressed in the form 


El-]tAi-ziy 

Z{zi,l) = l- 




(17.24) 


Comparing this with (17.23) we can make the identifications 


d-1 


d-1 


1 -^Sj{l - ZiY = ai = l - ^ 7 rf_j(l - ZiY, 
i=i i=i 

d-1 d-1 

'^tj{i- ziY = SjUj+i 
j=l 2=1 

d—1 d—1 j—1 

= 5; 5,(1 - z,y (17.25) 

j=l j=2 i=l 


which imply 




Sj = jd-j, j = 1,..., d-1, ti = di, tj = Sj-'^ 5ad+i-j j = 2 ,..., d-1. (17.26) 


2=1 


Given the coefficients s and t we can inductively uniquely determine the coefficients 7 
and 6 : 


j 

7i ^d—jj j 1,..., d 1 di ti, 5j tj ^ ^ 5iSj—i j 2,..., d 1. (17.27) 

2=1 

One can see from this analysis that there can be more than one pruned subspace 
collection associated with a rational function Z{zi, 1). It may happen that one pruned 
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Z{n) subspace collection gives rise to polynomials 1) = /(^i, l)r{zi, 1) and q{zi, 1) = 
g{zi, l)r{zi, 1) while another pruned Z(n) subspace collection gives rise to polynomials 
p'{zi,l) = f{zi,l)r'{zi,l) and q'{zi,l) = ^(zi, l)r'( 2 :i, 1 ), so that both give rise to the 
same function Z{zi,l). However there is a one-to-one correspondence if the pruned 
subspace collection is such that the polynomials p{zi,Z 2 ) and q{zi, Z 2 ) have no factor in 
common, and this correspondence is given by the above algorithm 


18 On the correspondence between certain rational 
functions of two variables and Z{3) subspace col¬ 
lections with m = 1 


In the case m = 1 and n = 3 can we uniquely recover a generic subspace collection (mod¬ 


ulo the linear transformations (13.9)) from knowledge of the rational function Z{zi,Z 2 , 1)? 


The answer is no, but let us hrst provide a counting argument which suggests that, at 
least in the generic case, we can recover the subspace collection up to a hnite number of 
possibilities. The counting argument is similar to that given in Section 29.2 of Milton 
( 2002 ) but here we do not assume that the subspaces are orthogonal. 

How many independent coefficients aa^a 2 a 3 are there in a polynomial 


p{Zi,Z 2 ,l)= 


a. 


^ 0,1 ^ 0,2 


aia2a3^1 ^2 i 


( 18 . 1 ) 


0-1,0.2,0.3 


that satishes 


“T U2 -r U3 




Without loss of generality, following Section 29.2 of Milton (2002), let us suppose that 
Pi > P 2 > Ps- With ai hxed in the regime 0 < oi < 1 -|- gi — p 2 , the constant 02 can take 
integer values from 02 = 1 + qi — ai — ps (where 03 = ps) to 02 = P 2 , that is, a total of 
P 2 +P 3 + cii — qi different values. With ai hxed in the regime p 2 < Oi < l + <?i—P 3 , 

the constant 02 can take integer values from a 2 = I + qi — ai — ps (where 03 = p^) to 
a 2 = I + qi — tti (where 03 = 0) that is, a total of ps -|- 1 different values. Finally, with ai 
hxed in the regime 1 -|- gi — p 3 < oi < pi, the constant 02 can take integer values from 
02 = 0 to 02 = 1 -l- gi — oi (where 03 = 0), that is, a total of 2 -|- gi — oi diherent values. 
Therefore the total number of coefficients in the polynomial is 


qi-P 2 <?i-P3 pi 

(P2 + P3 + ai - gi) + Y1 (^3 + 1) + Y2 (2 + gi - «i) 

ai=0 ai=l+iJi-P2 cil=l+<?i-P3 

= (gi -P 2 + i)(p 2 + P3 - gi) + ^(gi -P 2 + i)(gi -P 2 ) + (P 2 -P3)(P3 + 1 ) 

+(pi +P3 - gi)(2 + P3) - ^((pi +P3 - gi)(pi +P3 -gi + 1 ) 

= fci + l, (18.3) 

where 

h = [2(1 + gi)g 2 - P? - P 2 - Ps + ^]/2, (18.4) 
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in which h = pi + pi + and q 2 = h — 1 — qi. These coefficients are not all independent 
since, from (16.12) the aa^a 2 ai must sum to one. Subtracting this constraint gives ki 


independent coefficients. 
Similarly in a polynomial 


q{zi,Z2,l)= ^ 

^1162,63 


-61 ^P2-b2 
2 ; 


(18.5) 


that satishes 


bi + b2 + bs — 1 + q2, 0 < tti < Pi < 1 + q2, i — 1,2, 3, Pbib2b3 — I 5 (18.6) 

61,62,63 


there are a total of 


k 2 = |2(1 + 92)91 - p5 - p 1 - pi + A]/2 


(18.7) 


independent coefficients. Hence the total number of independent coefficients in the ra¬ 
tional function 

P{Z 1 ,Z 2 , 1 ) 


Z{zi,Z 2 , 1) = 


(l{Zi,Z 2 ,l) 


(18.8) 


IS 


ki + k 2 = (1 + qi)q 2 + (1 + q 2 )qi - pI - pI - pI +h = - pi - pi-pi-ql - ql- (18.9) 

Now how many parameters describe a Z{n) subspace collection, when the spaces U, 
£, J, Vi, V 2 , and "Ps have dimensions 1, qi, q 2 , pi, P 2 , and ps, with 1 + qi + q 2 = 
Pi + P2 + Ps = h? Let wi,W 2 ,... ,w/i be a basis for PL with wi in U, W 2 ,W 3 ,... ,Wgj+i in £, 
and Wg^+ 2 , Wqj+ 3 ... ,w/j in J'. Recall that it requires s{d — s) parameters to describe the 
orientation of a subspace of dimension s in a space of dimension d. Therefore, it requires 


Pi{h - Pi) + {h- P2)P2 + {h- P3)P3 = h'^ -pI-pI-pI (18.10) 

parameters to describe the orientation of the subspaces Vi, V 2 , and P 3 with respect to 
this basis. However some of these subspace collections are equivalent, linked through 
transformations of the form (13.9). If respect to this basis C is represented by a matrix 
with block form 

/c 0 0 \ 

C = 0 Cl 0 , (18.11) 

Vo 0 C 2 / 

where c is a scalar, while Ci and C 2 are qi x qi and q 2 x q 2 matrices, then it will leave 
the subspaces U, £ and J unchanged. The transformation C = al leaves all subspaces 
unchanged for any scalar a 7 ^ 0 , and so to factor out such trivial transformations we 
should choose c = 1. The number of remaining independent parameters in C is then 


ql + Subtracting these from (18.10) we see that the number of parameters describing 
the Z{n) subspace collection is 


1,2 2 2 2 2 2 7,; 

h Pi P 2 P 3 qi q2 — ki k2- 


(18.12) 


The precise agreement between the number of coefficients in the rational function and 
the number of parameters describing the Z{n) subspace collection is curious (since it 
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holds for all gi, q 2 , pi, P 2 , and ps, with 1 + gi + g2 = Pi + P 2 + P 3 = h). Despite this 
coincidence we now show that it is not possible to uniquely recover a generic subspace 


collection (modulo the linear transformations (13.9)) from knowledge of the associated 
rational function Z{zi, Z 2 , 1). 

Let us consider a subspace collection with h = 5, qi = q 2 = 2, pi = p 2 = I, ps = 3 
giving fci + /c2 = 6 according to the formula (18.9). Given u G W we choose as our basis 
the vectors 


Vo = U, Vi = AiU, V2 = A 2 U, V3 = riAiU, V4 = riA2U, 

with the closure relations 


(18.13) 


A1V3 = 7 iVi, A2V3 = 72 V 2 , AiV4=73Vi, A2 V4 = 74Vi, 

ToVi = ^iVo, roV2 = (52 Vo, (18.14) 

expressed in terms of the 6 parameters 7i, 72,73,74, (5i, and 62 which describe the sub¬ 
space collection. The question is: can one uniquely recover these six parameters from 
Z{zi, Z 2 , 1)3 Although the following analysis extends easily to the case of arbitrary 71 
and 74 let us assume, for simplicity, that 71 = 74 = 0 and ask whether one can recover 
the remaining four parameters. The held e -|- E must have the expansion 


e -|- E — Uo^o T ^1^3 -|- CI2V4, 


(18.15) 


and consequently, setting z^, = 1^ 

j + J = [I -|- (zi — l)Ai -|- (z2 — l)A2](e -|- E) 

= aovo aiV3 a2V4 + (zi - l)(ao a273)vi -h (z2 - l)(ao ai72)v2(18.16) 

Since ri(j + J) = 0 we arrive at the equations 

0 = aiV3 -h 02V4 (zi - l)(ao + a273)v3 + (z 2 - l)(ao + ai72)v4, (18.17) 

implying 

+ (d ~ l)(flo + ^273) = 0, 02-1- (z2 — l)(oo -l- O172) = 0. (18.18) 

These equations have as a solution, 

ao = 1 - (zi- l)(z2 - 1)7273, 

ai = 73(d - 1)(^2 - 1) - (zi - 1), 

02 = 72 (d - l)(z2 - 1) - (z2 - !)• (18.19) 

Since 


ro(e-|-E) — OqVq, ro(j-|-J) — [oo-|-(2:i —l)(oo-|-0273)^i-l-(2:2~l)(flo+®i72)^2]vo, (18.20) 


we obtain 


Z(Zi,Z 2 , 1) 


^ (zi - l)(oo 0273)(^1 + (Z2 - l)(oo + Oi72)<^2 

Oq 

^ ^l(zi - 1) - 73'5i(Zi - 1 )(Z 2 - 1) + '^2(^2 - 1) - 72^2(D - 1)(^2 " 1) 

1- (Zi - 1)(Z2 - 1)7273 


(18.21) 
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Given this function we can uniquely determine (5i and S 2 from the coefficients of {zi — 1) 
and {z 2 — 1) in the numerator. Also from the coefficients of (.^i — 1)(^2 ~ 1) in the 
numerator and denominator we can uniquely determine 


h — 7273 ) G — 735 i + 72 ^ 2 , 


(18.22) 


in terms of which there are two possible values of 72, namely 


72 = 


ts ± - 3ti6i52 

26 i 


(18.23) 


Thus we cannot uniquely recover the subspace collection parameters from Z{zi, Z 2 , 1). 


It remains an open question, raised at the end of Section 29.2 of Milton (2002), as 
to whether in general one can uniquely recover the subspace collection parameters when, 
with respect to some inner product, the subspaces U, £ and J' are mutually orthogonal, 
and the subspaces Vi, V 2 and V 3 are mutually orthogonal. These orthogonality con¬ 
straints overdetermine the system of equations needed to recover the subspace collection 
parameters which provides some hope that we can recover them. It would be useful if 
one could uniquely recover the subspace collection parameters (the weight and normal¬ 


ization matrices introduced in Milton, 1987a, 1987b) from say the effective conductivity 


(T*((Ti, (72, cTs) of an isotropic composite of three isotropic phases having conductivities 
(Ji, (J2, and (73 as then one could obtain the effective response tensor for coupled held 
problems. We will see in Chapter 9 of this book (Milton 2016) that the effective re¬ 


sponse tensor just depends on the weight and normalization matrices for the uncoupled 
conductivity problem. 


19 Visualizing the poles and zeros of functions asso¬ 
ciated with orthogonal Z{3) subspace collections 
when m = 1 


For scalar functions Z{zi, Z 2 , Z 3 ), associated with orthogonal Z{3) subspace collections, 
satisfying the homogeneity, Herglotz, and normalization properties, the trajectories of 
their poles and zeros in {zi,Z 2 ,Z 3 ) space, with zi, Z 2 , and z^ taking real values, have 
a beautiful visualization as trajectories on three interlinked hexagons: To obtain this 
visualization we follow Appendix C in Nicorovici, McPhedran, and Milton (1993): see 
also hgure 5 in that paper. 

First note that if we set Z 3 = 1, then the poles and zeros of Z{zi, Z 2 , 1) lie in one of 
the three quadrants: 


The quadrant Zi < 0, Z 2 > 0] 


• The quadrant -^2 < 0, zi > 0; 

• The quadrant Zi < 0, Z 2 < 0. 

Of course we can visualize the pole and zero trajectories by plotting them in this plane, 
but this has the disadvantage that the three variables Zi, Z 2 and Z 3 are not treated in 
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a symmetric way, and the disadvantage that its hard to see what is happening when zi 
and/or Z 2 is large, and it is hard to see what is happening near the origin zi = Z 2 = 0 
since the trajectories can bunch up there. To get around this we map each of the three 
quadrants to a hexagon. Given a quadrant, the point Zi = Z 2 = 0 gets blown up to form 
one edge of the hexagon; the two edges of the quadrant where zi or Z 2 is zero, but not 
the other, get mapped to two other edges of the hexagon; the two “boundaries” of the 
quadrant where \zi\ or \z 2 \ is inhnite but other is hnite get mapped to two more edges of 
the hexagon; hnally Zi = Z 2 = 00 gets mapped to the hnal sixth edge of the hexagon. We 
remark that just as a pole trajectory can cross from one quadrant to another, so too can 
it jump from the boundary of one hexagon to the corresponding point on the boundary 
of another hexagon. 

To be more precise, we introduce the three variables 


= 


1 + ^2/^3!’ 


to — 


1 


1 + ks/ Zl\ ’ 


t.s — 


1 + 1^1/^2! 


(19.1) 


Clearly (G, t 2 , t^) takes values in the unit cube. It is conhned to a surface within the unit 
cube as the three ratios 12:21/12:31, l^sl/l^il and |zi|/|.22| are not independent, but have 
product 1. The next step is to map these three variables onto three variables si, S 2 and 
S3 lying in the plane si + S2 + S3 = 0 using the projection 


Si — 2 ti — t2 — ts, S2 — 2^2 — G — S3 — 2^3 — ti — t2. 

Finally, we map these down to the x-y plane: 

y ^ (si + 2 s2)/v^. 


(19.2) 


(19.3) 


Some normalization is needed, so in the hexagon where zi is negative and Z 2 and Z 3 
are positive, we plot Z(zi, Z 2 , z-^j^fz^\ in the hexagon where Z 2 is negative and Zi and 
Z 3 are positive, we plot Z{zi, Z 2 , z^)/yjz^z^] and in the hexagon where Z 3 is negative and 
zi and Z 2 are positive, we plot Z[zi^ Z 2 , Z 3 )/^ziZ 2 - 

Figure louses this approach to visualize the pole trajectory of a function Z{zi,Z 2 , Z 3 ) 
associated with a Z(3)-subspace collection 


l~i — lA ® S ® iH — "Fi © 1^2 © ”^3, 


(19.4) 


where in this example Ti is 12-dimensional; U is one-dimensional; Vi is 3-dimensional; V 2 
is 6-dimensional; V 3 is 3-dimensional. Note that as the subspace collection does not need 
pruning, the dimensions of Vi, V 2 , and V 3 can be immediately read off from the hgure by 
simply counting the number of pole paths on each hexagon: hgures (a), (b), and (c) have 
3, 6 and 3 pole paths corresponding to the dimensions of Vi, V 2 , and V 3 , respectively. 
To understand this, hrst recognize that when Z 2 and Z 3 are hxed, and real and positive, 
Z{zi, Z 2 , Z 3 ) is a Herglotz function of Zi taking real positive values when zi > 0. Thus all 
its poles must be simple and located on the negative real Zi-axis, i.e. on the hexagon (a). 


Also because the subspace is pruned Ai(W © J') can be identihed with Vi (Section 16), 
and hence the matrix Ci representing Ai(ro + r2) has rank pi. Then as Ci and 
have equal rank (this well-known fact can easily be seen by showing that they have the 
same null-space), and as the subspace collection is orthogonal, it follows that the matrix 
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Bi = CfCi representing (ro+r 2 )Ai(ro + r 2 ) has exactly rank pi. Similarly the matrix 
Ai representing (ro + r 2 )Ai(ro + r 2 ) has exactly rank pi. Therefore the sum over ai in 
the numerator in (16.3), goes up to ai = pi, while the sum in the denominator in (16.11), 
goes from 0 up to bi = pi or (when all the coefficients zero) to = Pi — 1: it 

cannot go only up to 6 i = pi — 2 , since as a function of zi, Z{zi, Z 2 , Z'^/^fz 2 zi with hxed 
2:2 > 0 and hxed Z 3 > 0 can only have a simple pole at 2:1 = 00 . When the sum over bi 
goes up to bi = pi, there are clearly p poles of the function Z{zi, Z 2 , Zj)l^z^Zj, on the 
hexagon as Z\ varies with hxed 2:2 > 0 and hxed z^ > 0. When the sum over bi goes up 
to bi = Pi — 1 , there are still p poles of the function Z{zi, 2 : 2 , 2 : 3 )/^^ 2^3 on the hexagon 
as zi varies with hxed 2:2 > 0 and hxed 2:3 > 0 provided we count the pole at = 00 . 



Figure 8: The pole trajectory of the function Z(2;i, 2:2, 2:3) as visualized using the rep¬ 
resentation using three interlinked hexagons. The hexagon in (a) corresponds to real 
values of (2:1,2:2, 2:3) where 2:2 and 2:3 have the same sign, but Zi has the opposite sign. 
The hexagon in (b) corresponds to real values of (2:1, 2:2,2:3) where Zi and 2:3 have the 
same sign, but Z2 has the opposite sign. The hexagon in (c) corresponds to real values 
of (^1, Z2, Z3) where zi and 2^2 have the same sign, but 2:3 the opposite sign. By superim¬ 
posing all three pictures one obtains (d) where the pole trajectory is like that of a billard 
ball bouncing around a hexagonal table, following curved paths. The zero trajectory is 
similiar, but for clarity we chose not to include it. Note that the dimensions 3 , 6 and 3 
of the subspaces Vi, V2, and V3 can be immediately read off from the number of paths 
crossing the hexagons in (a), (b) and (c). 


The dimension q2 of the subspace can also generically be read off from the pole 
trajectories on the three hexagons. Consider the edge joining two of the hexagons that 
corresponds to the values 2:2 = 0 , and 2:3 = 1 with 2^1 < 0 varying. Then the only 
coefficients ^^at can contribute to the denominator in ( 16 . 11 ) are those with p2 = 
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b 2 . The first constraint in (16.10) then implies 


6l + 63 = 1 + ?2 - P2. 


(19.5) 


So bi can only range from 0 np to the maximnm of pi and l + q 2 —p 2 = Pi+P 3 — qi- Note 


that according to the ineqnality (15.6), qi > ps — 1 so 1 + ^2 — P 2 conld be as large as 
Pi + 1. If there are less than pi pole trajectories crossing this edge joining the hexagons, 
the nnmber of these crossing pole trajectories shonld generically allow ns to determine 
q 2 and hence gi, assnming pi, p 2 and pa have been determined from the nnmber of pole 
trajectories on each hexagon. If there are exactly pi pole trajectories crossing the edge 
then q 2 conld be pa or pa + 1. To determine which it is (or as an additional check on the 
valne of ^ 2 ) we conld look at pole trajectories, or zero trajectories, crossing other edges 
where the hexagons meet. 

This visnalization may be nsefnl in hnding other topological featnres of the trajec¬ 
tories, which hopefnlly conld be connected with topological featnres of the snbspace 
collections. 


20 Normalization operations on subspace collections 


Rational fnnctions of a single variable may be expanded in continned fractions, which 
incorporate snccessively higher and higher order terms in the series expansion of the 
fnnction abont a point. The analogons procednre with snbspace collections is achieved 
throngh normalization and rednction operations, snbject to some technical assnmp- 
tions. The associated fnnctions are then linked, and provided the technical assump¬ 
tions hold at each level, these links provide continued fractions for multivariate functions 
Z(zi, ^ 2 , • • •, Zn) and Y(zi, ^ 2 , • • •, Zn) incorporating matrices of increasingly high dimen¬ 
sion at each level in the continued fraction. 

The normalization and reduction operations are discussed in this and the next section. 
For more insight, in the case where the subspaces in the direct sums are orthogonal (see 


Milton 1987a, 1987b and Sections 19.2, 20.6 and 29.5 in Milton 2002). 


Normalization reverses extension. Given a subspace collection 


jC = £' Q)J' = V®Vi®V2 




( 20 . 1 ) 


dehne 


n = Vi®V 2 ®---®Vn, s = 8'nn, j = j'nn, 
u = n2r;v = n2(i - r'2)v = n2r'2V, 8 = r;v, j = r'2V, ( 20 . 2 ) 

where and r2 are the projections onto 8 ' and J\ and 112 is the projection onto H. 

We assume that the Y-problem has a unique solution when L = I for Ji G V given 
El G V. In other words, we assume that the equations 

El -|- E 2 G 8 \ Ji -|- J 2 G J 2 = E 2 , El, Ji G V, E 2 , J 2 £ 77, 

El + E 2 G 8 ', Ji + J 2 e J7', J 2 = E 2 , Ii G V, E 2 , J 2 G 77, (20.3) 

imply Ji = Ji. Subtracting these equations we see that 

E = E2 — E2 G 8 ' ^ J = Ji-|-J 2 — Ji — J 2 ^ 1 J2 — J2 = E. (20.4) 
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These imply 


E G 'H, E = J — V, where v = Ji — J^. (20.5) 

The uniqueness assumption means that these equations imply v = 0 (and if v = 0 then 
necessarily E = J = 0 since E' and J' have no vector in common). The relation E = J —v 
with E G n "H implies 

E = -r'lV, (20.6) 

which will only have the trivial solution v = 0 if and only if 

Hn£ = 0 and Vnj' = 0, (20.7) 

where the latter guarantees that r)v = 0 implies v = 0. 

We also assume that the E-problem has a unique solution when L = I for Ei G V 
given Ji G V. By similar analysis this is satished if and only if 

Hnj = 0 and VnT' = 0. (20.8) 

We now establish that 

= E®J = V ®U. (20.9) 

First note that V and U. have no vector in common since U G H, and similarly S and 77 
have no vector in common since £' Hj'' = 0. Clearly W contains V. To show it contains 
U notice that 


u = U2 T[v = (I - ni)r;v c r;v © nir;v c ^ © v c w. 

Together these imply V © 7/ C W. Finally we have 

s = T[v = (Hi + n2)r;v c nir;v © n2r;v c v ©w, 


( 20 . 10 ) 


( 20 . 11 ) 


and similarly 77 C V ©7/. Together these imply W C V ©7/, establishing (20.9). 

If V has dimension m then S must also have dimension m since otherwise r)v = 0 
for some nonzero v G V, implying v = r 2 ’'^ which on ly ha s the solution v = 0 since 
V n 77 ^ = 0. Similarly 77 must have dimension m and (20.9) then implies U must have 

W=U®S, ( 20 . 12 ) 


dimension m. The hrst condition in (20.7) implies 


since lA Gl-i and £ have no vector in common and are m-dimensional spaces contained 
in the 2m-dimensional space W. Now any vector E' G £' has the unique decomposition 


E' = e; + p, e; g V, 

P G 77, 

(20.13) 

and according to (20.12) E) has the unique decomposition 


E]^ — —e + E, e G 77, 

E G 

(20.14) 

So we have the decomposition 

E' = E + E, 


(20.15) 
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where 


E = P- e = 


(20.16) 


Also the hrst condition in (20.7) implies £ and £ G H have no vector in common, so the 

(20.17) 


decomposition is unique. Therefore we conclude that 

£' = £®£, 


and similarly the hrst condition in (20.8) implies 

J' = J®J. 


These and (20.9) imply 

K. = V®'H = £®£®J®J = V®U®£®J, 


(20.18) 


(20.19) 


and since U, £ and J are all contained in 77 we conclude that 


n = u®£®j = Vi®V2 


V 

f n. • 


( 20 . 20 ) 


Now a given E( G V has the unique decomposition (20.14). This dehnes the nonsin¬ 
gular operator K. : V U such that e = KE). (It is nonsingular because V and £ G £' 
have no nonzero vector in common.) Now given e, consider the solution to 


n 

e, j G 77, E G T, J G J7, j + J = L(e-|-E), where L = ^ (20.21) 

i=l 


where A* is t he pro jection onto Vi, and from the dehnition of Z, j = Ze. Since the second 
condition in (20.7) implies V and J have no vector in common we have 


W = V © 77, 


( 20 . 22 ) 


and consequently any j G 77 has the decomposition 

j = —Ji + J, J) G V, J G 77, (20.23) 

which dehnes the nonsingular operator M ; 77 —)■ V such that J) = Mj. Dehning 



E 2 — e + E, J 2 — j + J, 

(20.24) 

we have 

E'^ + E 2 = E) + e + E = E + E G £' , 
j; + J'2 = j;+j + J = J + JG77', 

(20.25) 

and 

j; = Mj = MZe = mzke;, 

which by dehnition of the associated E-function implies 

(20.26) 


Y(zi, Z2, • • • , Zn) = MZ(zi, Zn)K. 

(20.27) 
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This is analogous to the relation (20.29) in [Milton (2002 ) obtained in the case where the 
subspaces are mutually orthogonal. 

In particular by letting zi = Z 2 = ■ ■ ■ = Zn = 1 we obtain 


Y(1,1,...,1) = MK. 


(20.28) 


If vi, V 2 ,. .. Vm are a basis for V, and we choose Kvi, Kv 2 ,. .. Kv^ as our basis for W 
then with these bases K is represented by the identity matrix K = I and (20.27) and 

Y (^ 1 , ^ 2 , • • •, ^n) = Y(l, 1,..., l)Z(zi, Z 2 ,..., Zn). (20.29) 


([20^ imply 


21 Reduction operations on subspace collections 

Extension is one way to go from a Z (n) subspace collection to a Y {n) subspace col¬ 
lection. Another way is through reduction, which has some features in common with 
normalization. Given a Z{n) subspace collection 

n = u®s®j = Vi®V2®---®Vn, ( 21 . 1 ) 

let To be the projection onto W, and let Aj be the projection onto Vj. Dehne 

1C = £ ® J, Vj = Vj n 1C for j = 1, 2,..., n, 

V = (I-ro)[AiW©A 2 W©---©A,W] c/C, Vj=AjU. (21.2) 

We now establish that 


W = © 7^2 © ■ ■ ■ © 'Pn — 77 © V. 


(21.3) 


First note that V and U have no vector in common since V G 1C, and similarly the 
subspaces Vj have no vector in common since Vj C Vj. Clearly W contains U since the 
projections Aj sum to the identity. To show it contains V note that 

V C AiU © A 2 U © • • • © AnU + ToiAiW © A 2 U © • • • © AnU] C W + W = >V. (21.4) 


Therefore we have that 7/ © V C W. The converse inclusion that W C 77 © V follows 
from the inclusion 

Vj = [To + (I - To)]AjU C 77 © V, (21.5) 


which establishes (21.3). Next, to establish that for all j, 

Vj = Vj®V'j, 


we need to assume that for all j 
and that 


VjnlC = 0, 


A,u = 0, u G 77 


only has the trivial solution u = 0, i.e. 

77 n {Vi © V 2 © ... © Vj—i © Tj+1 ©... © Vn) = 0. 


( 21 . 6 ) 

(21.7) 

( 21 . 8 ) 

(21.9) 
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These conditions imply that 

and hence that 

which in turn implies that 


U = ToA^W, 

U C AjU © (I - To)AjU, 
U c Vi + V. 


Then any vector P G Vj has the unique decomposition 

P = u + K, with u K G /C, 


and according to (21.12), u has the unique decomposition 

u = V + P with V G V, P G Vj, 


( 21 . 10 ) 

( 21 . 11 ) 

( 21 . 12 ) 

(21.13) 

(21.14) 


which is unique because V C /C and Vj have no nonzero vector in common. Therefore P 
has the unique decomposition ^ 

P = P + P', (21.15) 


(21.16) 


where 

P' = V + K = P - P G Pj n /C = Pj. 

This decomposition and the fact that (21.7) implies Vj and P) C /C have no vector in 
common establishes (21.6). 

So we deduce that 

'H=W©^©d7 = Pi©p2©---©P„©P(©P2©---©P(, 

= W © V © P( © p2 © • • • © P(„ 

and since the P), j = 1, 2,..., n are all contained in /C it follows that 

/C = ^ © JT = V © P( © p2 © • • • © P(,. 

Now suppose that given e G W we can solve the equations 

j + Ji = L(e + El), Ji = —YEi, e, j G W, Ei,JiGV, 


(21.17) 


(21.18) 


(21.19) 


where Y is the Y-operator associated with the subspace collection (21.18). From the 
Y-problem we have 


where 

Since 


E — El + E 2 g£^, J — J1 + J 2 G 77 J2 — LE2, E 2 , J2 G l-i\ 

P' = P( © P2 © • • • © P(,. 


j + Ji + J2 — L(e + El + E2), 
we see that these helds solve the Z-problem 

e, j G 7 /, E G T, J G 77, j + J = L(e © E), 


( 21 . 20 ) 

( 21 . 21 ) 

( 21 . 22 ) 

(21.23) 
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and by definition j = Ze. To solve (21.19) let Hi be the projection onto V. Then ( |21.19[ ) 
implies 

-YEi = niL(e + niEi), 


giving 

El = -ni(Y + niLHij-^niLe, 

where the inverse is to be taken on the subspace V. It follows that 

j + Ji = Le - Lni(Y + niLHi) ^HiLe, 

implying 

Z = EoLEo - EoLni(Y + niLni)-^niLEo. 

This formula is analogous to that given in (29.12) of Milton (2002| ). 

To obtain a more explic it wa y of writing (21.27) let us suppose we are 
ui, U 2 ,..., Urn of U. Since (21.8) only has the trivial solution u = 0 each 


dimension m. It then follows from (21.3) that V has dimension m{n 
i = 1,2,..., n — 1, (21.3) implies AjUj has the unique decomposition 


(21.24) 

(21.25) 

(21.26) 

(21.27) 

given a basis 
space Vj has 
1 ). Also, for 


AiUj = ^ Wp-fcUfc + Vij, G V, 


(21.28) 


for some set of constants tcp-fc. To show that the vectors Vjj, which number m{n — 1), are 
independent, let us suppose 


n—1 m n—1 m m 

WijkUk). (21.29) 

2 = 1 j=l 2=1 j = l k = l 


By letting A„ act on this equation and taking into account that (21.8) only has the trivial 
solution u = 0 we see that 


n—1 m m 

EEE ^ij'^ijk^k 0. 

2=1 j=l k=l 


(21.30) 


Then substituting this in (21.29) and letting Aj, i ^ n, act on (21.29) and again taking 
into account that (21.8) only has the trivial solution u = 0 we obtain 


Ci3^3 = 
i=i 


(21.31) 


which shows that all the must be zero. Therefore let us take the vectors v^- as our 
basis for V. 

The identities 


Il-iAjEoUj 


rcAjEgUj ^ ^ Wjjk^ki 


k 


(21.32) 
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which follow from (21.28) then gives the matrix representations for niAiFo and roAiFo 
in these bases, when i ^ m. Using the fact that A„ = I — obtain 


n—1 n—1 

ToLFo = ^nro + — 2;n)roAjFo, FIiLFo = — ^n)niAjFo. (21.33) 

i=l i=l 

Now for p ^ n (and i ^ n) ( |21.28[ ) implies (no sum over p) 

k 

^ ^ijk')i^pk “1“ ^ ^ ^pkq^q) ■ (21.34) 

k q 

Thus we deduce 

roApFIiVjj ^ ^j i^^pi^kj '^ijk') ^ ^ Wpl^qVlq^ 
k q 

niApFIj^Vjj ^ ^ j i^^pi^kj ^ijk^^pkj (21.35) 

k 

which gives the matrix representation for the operators FoApFIi and FIiApFIi in these 
bases {p ^ n), in terms of which we obtain the representation for the operators 


n—1 


n—1 


FoLFIi — — ^ji)FoApFIi, FIiLFIi — ^„FIi + ^^( 2 ;^ — ^„)FIiApFIi. (21.36) 

P=1 


Thus all the matrices representing the operators entering (21.27), aside from Y, only 
depend on the parameters Wijk and these parameters can be obtained from the represen¬ 
tation in the basis ui, U 2 ,..., of Z when the differences z* — i = 1, 2,..., n — 1 are 
small. To hrst order in these differences, (21.27),(21.33), and (21.36) imply 


n—1 


ZUj Ri ZnUj + '^{Zi - Zn)'^WijkUk. 


(21.37) 


2 = 1 


Thus knowing this expansion one can recover all the parameters Wijk- 


22 “Continued fraction expansions” of subspace col¬ 
lections. 

The idea to developing the continued fraction is that by a succession of reduction and 
normalization operations one obtains a series of recursion relations 

Z = FoLFo - FoLni(Y + niLni)-^niLFo, (22.38) 

Y = (22.39) 

z(^) = (22.40) 

Y(i) = (22.41) 

z(2) = - FE,^^L(2)nS^^(Y(2) + (22.42) 
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and so forth, until the dimension of the remaining space goes to zero, or until one (or 
more) of the assumptions necessary to proceed with the normalization or reduction op¬ 
eration does not hold. By substituting (22.39) in (22.38), then substituting (22.40) in 
the resulting expression, and subsequently substituting (22.41) in this expression, and so 
on, one develops the continued fraction expansion for Z incorporating the variables zi, 
Z 2 , ..., Zn and, as one goes down the continued fraction, information contained in the 
series expansion (15.1) at successively higher and higher levels of truncation. We do not 
address in this book whether one can go ahead with the continued fraction expansion 
(and if so how) when the assumptions made to proceed with the normalization or reduc¬ 
tion operation do not hold. In the process of developing the continued fraction through 
reduction and normalization operations, one could at those steps where one is dealing 
with a y-problem make any desired reference transformation as described in Section 12. 
In this way one incorporates information at the subspace collection level that corresponds 
at the function level to known values of the function, and derivatives, at various points. 

Such continued fraction expansions form the basis of the held equation recursion 
method for bounding the effective moduli of composites (Milton and Golden 1985 Milton 


1987a, 

1987b 

1991 

Clark and Milton 1994 

Clark 1997 

and Chapter 29 of ( 

Milton 2002 


in the abstract theory of composites as described in Chapter 2 of this book (Milton 2016): 
see also Section 9.10 and Chapter 10 of Milton (2016)). The basic idea, at least when we 
have an orthogonal subspace collection, is that crude estimates or bounds on the operator 
ZO) or at some intermediate level j give through the above recursion relations good 
approximations or tight bounds on Z or Y incorporating the parameters that enter the 
recursion relations at the different levels up to level j (obtained from series expansions 
up to a given order of the solutions of the Z-problem or Y-problem). 
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